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3; Facility: 
VAX=-11 Instruction Emulator 

; Abstract: 
The routines in this module emulate the VAX-11 packed decimal 
instructions that perform arithmetic operations. These procedures can 


be a part of an emulator pogeege or can be called directly after the 
input parameters have been loaded into the architectural registers. 


The input parameters to these routines are the registers that 
contain the intermediate instruction state. 
| 


These routines run at any access mode, at any IPL, and are AST 
reentrant. 


; Author: 
Lawrence J. Kenah 
; Creation Date 

; 19 October 1983 | 
: Modified by: | 
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v01=003 LJK0037 Lawrence J. Kenah 17-Jul-1984 
Fix two minor bugs in exception handling code ehet caused 
MULP and DIVP tests to generate spurious access violations. 


v01-002 LJK0024 Lawrence J. Kenah 21-Feb-1984 
Add soee to handle access yiolations. Perform minor cleanup. 
Eliminate double use of R10 in MULP and DIVP. 


v01-001 LJK0008 Lawrence J. Kenah 19-0¢ t-1983 
The emulation code for ADDP4, ADDP6, SUBP4, SUBP6, MULP and 
DIVP was moved into a separate module. 
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Declarations P-1984 34 CEMULAT.SRCJVAXARITH.MAR; 1 

. SUBTITLE Declarations 
; Include files: 

-NOCROSS 3; No cross reference for these 
» ENABLE SUPPRESSION ; No symbol table entries either 
ADDP4_DEF ; Bit fields in ADDP4 registers 
ADDP6_DEF ; Bit fields in ADDP6 registers 
DIVP_BEF ; Bit fields in DIVP registers 
MULP_DEF ; Bit fields in MULP registers 
SUBP4_DEF ; Bit fields in SUBP4 registers 
SUBP6_DEF ; Bit fields in SUBP6 registers 
SPSLDEF ; Define bit fields in PSL 
SSRMDEF ; Define arithmetic trap codes 
-DISABLE SUPPRESSION ; Turn on symbol table again 
CROSS ; Cross reference is OK now 


; Symbol definitions 


The architecture requires that R4 be zero on completion of an ADDP6 or 
SUBP6 instruction. If we did not have to worry about resterstng 
instructions after an access violation, we could giapty zero the saved 
R4 value on the code pecs that these two instructions have in common 
before they merge with the ADDP4 and SUBP4 routines. The ability to 
restart requires that we keep the original R4 around at least until no 


SOOCOCCOCOCOCSCOSOSCOSCOCSOSOSOSOOSCSOSCOSOCOOOSOSOOOSOOOOOOOOoOoOO 
SOoOOCSCCCOOSOOCOOSOOSOOOOOOSOSOOSOOOCOOOSOOOOOOOSoOO 
COOOCCOCOCCOCOCOCOSOSGSOCOCCCOOCCOOCO OO OCC OO OOOO OOOCOOOOOOoOO 


OOOCCOCCSOOSOOOOCOOOCOOOOCOOOOOCOOOOOOOOoOO 
RPINIMNININININS 2 9 9 OF HQOOODOOOOOOOOOOO0 OOOO OW 00009 OOO NINN NNN 
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1 more access violations are possible. To accomplish this, we store the 
1 fact that R4 must be cleared on exit in R11, which also contains the 
1 evolving condition codes. We use bit 31, the compatibility mode bit 
! because it is nearly impossible to enter the emulator with CM set. 
0000001F ADD_SUB_V_ZERO_R4 = PSLSV_CM 
: : External declarations 
: DISABLE GLOBAL 
1 -EXTERNAL - 
1 DECIMALSBOUNDS_CHECK,- 
1 DECIMALSBINARY~TO_PACKED_TABLE,- 
0 1 DECIMALSPACKED_TO_ BINARY TABLE,- 
888 1 DECIMAL SSTRIF ZEROS RO. ” 
805 : DECIMALSSTRIP_ZEROS_R2_R3 
0000 1 EXTERNAL - 
0000. = VAXSDECIMAL_EXIT,- 
0 6 1 VAXSDEC IMAL -ACCVIO,- 
1 VAXSREFLECT_TRAP,- 
0 1 VAXSROPRAND 
48 1 : 
3 0 : 3; PSECT Declarations: 
8 ' . DEFAULT DISPLACEMENT , WORD 
0000 $0 1 -PSECT _VAXSCODE PIC, USR, CON, REL, LCL. SHR, EXE, RD, NOWRT, LONG 


dl 
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8000 189 BEGIN_MARK_POINT 
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: - SUBTITLE VAXS$SUBP6 = Subtract Packed (6 Operand Format) 
; Functional Description: 


In 6 operand format, the subtrahend string specified by the subtrahend 
length and subtrahend address operands is subtracted from the minuend 
str ng specified by the minuend length and minuend address operands. 
The difference string specified by the difference length and difference 
address operands is replaced by the result. 


; Input Parameters: 


sublen.rw Number of digits in subtrahend string 
subaddr.ab Address of subtrahend string 
minlen.rw Number of digits in minuend string 
minaddr.ab Address of minuend string 

diflen.rw Number of digits in difference string 
difaddr.ab Address of difference string 


: Output Parameters: 

; = 0 

: Ri = Address of the byte containing the most significant digit of 
; the subtrahend string 

: R3 = Address of the byte containing the most significant digit of 
: the minuend string 


Address of the byte containing the most significant digit of 
the string containing the difference 


; Condition Codes: 
N <- difference string LSS 0 


Z <- difference string EQL 0 
V <- aaa over flow 


Register Usage: 


COOOOCCOSGOCOSOSOOCOSCCOCSOOOOSOOOCOOOOOoOoO: 


This routine uses all of the general registers. The condition codes 
are recorded in R11 as the routine executes. 


- ENABLE LOCAL _BLOCK 


VAXSSUBP6: : 
PUSHR #*M<RO,R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,R11> ; Save the lot 
MOVZBL #1,R9 ; Indicate that this is subtraction 
BRB 10$ : Merge with ADDP6 code 
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. SUBTITLE VAXSADDP6 = Add Packed (6 Operand Format) 
Functional Description: 


In 6 operand format, 


he addend 1 string specified by 
Length and 


adden ddress operands is added to the addend 
specified by the addend Length and addend 2 address operands. 
string specified by the sum 
by the result. 


; Input Parameters: 


RO = addilen.rw Number of digits in first addend string 
R1 = addiaddr.ab Address of first addend strin 

R2 = add2len.rw Number of digits in second addend string 
R35 = addéaddr.ab Address of second addend string 

R4 = sumlen.rw Number of digits in sum string 

R5 = sumaddr.ab Address of sum string 


; Output Parameters: 


RO = 0 


@ 
par) 
" 


Address of the byte containing the most significant digit of 
the first addend string 


Address of the byte containing the most significant digit of 
the second addend string 


= 
Address of the byte containing the most significant digit of 
the string containing the sum 


~ 
uw 
" 


; Condition Codes: 


N <= sum string LSS 0 
Z <= sum string EQ. 0 
V <- lente overflow 


; Register Usage: 


This routine uses all of the general registers. The condition codes 


are recorded in R11 as the routine executes. 


VAXSADDP6: : 

PUSHR #*M<RO,R1,R2,R3,R4,R5,R6,R7,R8B,R9,R10,R11> 
CLRL : This is addition 
oor RG : Insure that R4 is LEQU 31 
MOVPSL R1 ; Get initial PSL 

3; Indicate that the saved R4 must be cleared on the exit path 


BBCS #ADD_SUB_V_ZERO_R4,R11,308 
BRB 30$ 


10$: 


the addend 1 


ength and sum address operands is replaced 


; Save the lot 


; Set bit and join common code 
:; In case we drop through BBCS 
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» SUBTITLE 


In 4 operand format, 


string spe 
operands an 


Parameters: 
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VAX$SUBP4 = Subtract Packed (4 Operand Format) 


Functional Description: 


cified by 
d the differe 


RO - sublen.rw 
R1 = subaddr.ab 
Re - diflen.rw 
R35 = difaddr.ab 


RO = 0 


Output Parameters: 


the 


the subtrahend string specified by suotrahend 
length and subtrahend address operands is subtracted from the 
difference length and difference address 


difference 


ce string is replaced by the result. 


A 
Number of digits in difference string 
Address of difference decimal string 


R1 = Address of the byte containing the most significant digit of 
the subtrahend string 


R2 = 0 
Rg = Address of the byte containing the most significant digit of 
the string containing the difference 


; Condition Codes: 


65 

66 N <= difference string LSS 0 

67 Z <- difference string EQL 0 

68 V <= decimal overflow 

$° C <- 0 

4 ; Register Usage: 

This routine uses all of the general registers. The condition codes 

2 Me are recorded in R11 as the routine executes. 

76 

77 VAXS$SUBP4: : 
OFFF 8F BB é 78 PUSHR #*M<RO,R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,R11> ; Save the lot 
59 01 QA 26 79 MOVZBL #1,R9 : Indicate that this is subtraction 

06 11 0029 280 BRB 20$ : Merge with ADDP4 code 
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SUBTITLE VAXSADDP4 ~ Add Packed (4 Operand Format) 
Functional Description: 
In 4 operand format, the addend string specified by the addend length 
and addend address operands is added to the sum string specified by the 
= Length and sum address operands and the sum string is replaced by 
ult. 


; Input Parameters: 


Rg - addlen.rw Number of gh HG. in addend string 
R1 = addaddr.ab Address of addend decimal string 
- sumaddr.ab Address of oun decimal string 


; Output Parameters: 


RO = 0 
R1 = Address of the byte containing the most significant digit of 
pe addend string 
= 


R35 = Address of the byte containing the most significant digit of 
the string containing the sum 


; Condition Codes: 


N <= sum string LSS 0 
i, <= sum string EQL 0 
4 <= decimal overflow 

{ <= 0 

; Register Usage: 

This routine uses all of the general registers. The condition codes 
are recorded in R11 as the routine executes. 


VAXSADDP4: : 


beDeste Testo ste Desde de Teste TesTesTes[esTesTesTesTesTestesTestecte Te sTes[esto[esTe[estesTesTosTes[esTecTesTosles) 
FW OOONOUES WN (OOD NAUE WIN 9 ODNAUES WN (OOD NA NEW O ODNOAUE UTP 


AINININIIPIPIPIPDPIPININININS 2 2 OO MH OOOOOOCOCOCOCOCOWOOOOWOOONO 


SSSssssssssssssssssssssssssssssssss 
SOOCOSOCOCOCOOCO COO SOSOC OOOO SOSO OOOO OOOO OOCOOOOOSOOOOOO 


OF FF ef BB PUSHR #*M<RO,R1,R2,R3,R4,R5,R6,R7,R8,R9,R10 R11> ; Save the lot 
9 D4 f CLRL = RO : This is addition 
1 3; The output string, described by R4 and RS, will be the same as the input 
1 3; string for ADDP4 and SUBP4. It is necessary to explicitly clear R4<31:16> 
! ; along this code path so MOVO R2,R4 will not always work. 
54 2g 3¢ 1 208: MOVZWL R2,R4 ; Set output size equal to input size 
a d0 4 MOVL R3,R5 3.2. and ditto for string addresses 
5B oO z MOVPSL Rif > Get initial PSL 
, : Indicate that the saved R4 will be restored on the common exit path 
00 58 #«1F~sCESS 33 BBCC #ADD_SUB_V_ZERO_R4,R11,30$ ; Clear bit and join common code 


| 
| 
Re - sumlen.rw Number of digits in sum string 
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§ - SUBTITLE ADDPx/SUBPx Common Initialization Code 


ALL four routines converge of this point and execute common initialization | 
code until a later decision is made to do addition or subtraction. 


R4 = Number of digits in destination string 
RS = Address of destination string 


R9 = Indicates whether operation is addition or subtraction 
0 => addition 
=> subtraction 


Www 


ROOD NAN ES WN @§ OW NAUNE WIN OS OD NAW EWN OS OD NAME WO OONOAUE Win 0 


R11<31> = Indicates whether this is a 4-operand or 6-operand instruction 
0 => 4-operand (restore saved R4 on exit) 
1 => 6-operand (set R4 to zero on exit) 


Sete Ge Se Ge Ge Ge Ge Ge Se Se Ge Se Se Ge 
* 


5B 046 00 04 F0 30$: NSV #PSL$M_7,#0,44,R11 ; Set Z-bit, clear the rest in saved PSW 
STABLISH_HANDCER ad 3; Store address of access 


ARITH_ACCVIO 3 violation handler 


ROPRAND CHECK R2 : Insure that R2 is LEQU 31 
MARK_POINT ADD_SUB_BSBW g 
BSBW DECIMALSSTRIP_ZEROS_R2_R3 ; Strip high order zeros from R2/R3 


ROPRAND_CHECK RO ; Insure that RO is LEQU 31 
MARK_POINT ADD_SUB_BSBW_0 
BSBW DECIMALSSTRIP_ZEROS_RO_R1 ; Strip high order zeros from RO/R1 


; Rather than es ges gentuse the already complicated logic dealing with 

; different length strings in the add or subtract loop, we will put the 
result into an intermediate buffer on the stack. This buffer will be long 

enough to handle the worst case so that the addition loop need only concern 

itself with the Lengths of the two input loops. The required Lengt is 1 

bytes, to handle an addition with a carry out of the most significant byte. 

We will allocate 20 bytes to maintain whatever alignment the stack has. 


CLRQ -(SP) 
CLRQ -(SP) 


CLRL -(SP) 

EXTZ #1,44,R4,R8 

ADDL R5 (SP) 

MOVAB 24(SPS,R 
The number of minus signs will determine whether the real operation that we 
perform is addition or subtraction. That is, two plus signs or two minus 
signs will both result in addition, while a plus sign and a minus sign will 


result in subtraction. The addition and subtraction routines have their own 
methods for determining the correct sign of the result. 


I 
E 
FFAE’ 30 


FFA3" 30 


3; Set aside space for output string 
; Worst case string needs 16 bytes 
3; Add slack for a 


PAREAANAAAANAAAAAAATE PERE 


Get byte count of destination string 
; Save high address end of destination 
Point R5 one byte beyond buffer 


OOOOCCOCOOOOoOoOoOooooooooo 


COOCOCSCSOCSOOGOOCOOOSOSSSOSOSCSOCSOSOSOSOSOOOOoOoOO 


~ 

m 
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For the purpose of counting minus signs. we treat subtraction as the 
addition of the negative of the input operand. That is, subtraction of a 
positive quantity causes bi sign to be remembered as minus and counted as 
a@ minus s an while subtraction of a minus quantity stores a plus sign and 
counts nothing. 


4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
5 
5 
5 
5 
P] 
5 
5 
5 
5 
5 
6 
6 
6 
6 
6 
6 
6 
6 
? 
? 
? 
? 
? 
? 
? 
? 
9 
9 


On input to this code sequence, RI distinguished addition from subtraction. 


kK 9 
VAX IMAL _ARITHMETIC = VAX=11 Packed Decimal Arithmetic Instr 16-SEP-1984 01:33:44 YAX/VMS Macro V04-00 P 1 
vourb00 itializati SEp-198¢ Ooaerse age 8) 


ADDPx/SUBPx Common Initialization Code 5-SEP-1 EMULAT.SRCJVAXARITH.MAR; 1 
7 3; On output, it conteing either 0, 1, or 2, indicating the total number of 
f 4 ; minus signs, real or implied, that we counted. 
56 50 o% 9} EF 7 9 EXTZV #1,#4,R0,R6 ; Get byte count for first input string | 
1 6 CO 7 9 ADDL R6,R1 ; Point R1 to byte containing sign 
56 61 FO 8F ; ; BICBS  wAB111 100007 RTT -RG R6 tains the si digit | 
s ° ° 3 contains the sign "' - 
8 4 8 07 : BLBS R9,35$ 3; Use second CASE 1} "subtrection 
0 Z § ; This case statement is used for addition 
404 CASE R6,TYPE=B,LIMIT=#10,<- ; Dispatch on sign digit 
6 405 G4! 3; 10 3 sign is aes ° 
0 4 $ 40$,- 3 11 => sign is ‘= 
0 4 50$,- ; \¢ => sign is ‘'+"’ 
080 $88 40$,- 3; 13 => sign is ‘'="' 
0 409 50$,- 3 14 => sign is ‘‘+" 
0 410 50$,- 3; 15 => sign is ‘'+" 
+ 411 > 
090 216 
494 ai7 3; This case statement is used for subtraction 
0090 415 35$: CASE R6, TYPE=B,LIMIT=#10,<- ; Dispatch on sign digit 
0090 $18 40$,- ; 10 => treat sign as ‘'-"' 
0090 41 50$,- ; 11 => treat sign as ‘'+" 
009 418 40$,- 3 \¢ => treat sign as ‘= 
09 419 50$,- : 135 => treat sign as ‘'+"’ 
090 420 40$,- : 14 => treat sign as ‘'="' 
0090 421 40$,- ; 15 => treat sign as ‘'~"' 
0090 4 ¢ = 
OOA 4 
59 01 DO O0A 424 40$: MOVL #1,R9 ; Count a minus sign | 
56 OD 9A OQOA 425 MOVZBL #13,R6 ; The preferred minus sign is 13 
05 11 pas : § BRB 60$ ; Now check second input sign 
59 D4 QOAB 428 50S: CLRL R9 ; No real minus signs so far 
56 OC 9A ~~ ? 4 MOVZBL #12,R6 ; The preferred minus sign is 12 
57 52 046 01 EF QOAD 431 60$: EXTZV #1,#4,R2,R7 ; Get byte count for second input string 
ss. rT. @ bes 4 § ADOL R7,R3 ; Point R3 to byte containing sign 
085 4 MARK_POINT ADD_SUB_ 24 ; 
57 63 FO 8F 8B ¥4 ? : BICB3 #*B111100007(R3),R7 3; R7 contains the sign ‘‘digit’’ 
GOA 4 § CASE R7,TYPE=B,LIMIT=#10,<- ; Dispatch on sign digit 
A 4 0$,- : 10 => sign is “+ 
OBA 438 0$.- : 11 => sign is "=" | 
A 439 to ; \ => sign is 
A 440 0$,- ; 13 => sign is ‘'-"" 
A 441 80$,- ; 14 => sign is ‘'+"' 
A rei 80$,- ; 15 => sign is *'+ 
sa 44 > 
CA 444 
59 D6 OOCA 445 708: INCL R9 ; Remember that sign wes minus | 
57 3 SA CC 44 MOVZBL #13,R7 ; The preferred minus sign is 13 
0 11 4 2 BRB 90$ ; Now check second input sign 
57 OC 9A D1 ra: 80$: MOVZBL #12,R7 ; The preferred minus sign is 12 
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2 90$: BLBC R9,ADD_PACKED 

2 i BRW SUBTRACT_PACKED 

455 - DISABLE LOCAL BLOCK 
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3; Even parity indicates addition 
; Odd parity calls for subtraction 


ae. 7 


ODA 4 Note that the correct choice of sign in all four cases is_the sign 
oe of the second input string, the one described by R2 and R3. 

OODA Input Parameters: 

QODA 


RO<4:0> - Number of digits in first input decimal string 
R1 - Address of least significant digit of first input 
decimal string (the byte containing the sign) 


~ rele - Number of digits in second input decimal string 
R Address of least significant digit of second input 
decimal string (the byte containing the sign) 


R4<4:0> = Number of digits in output decimal string . 

R5 - Address of one byte beyond least Significant digit of 
intermediate string stored on the stack 

RE<d:> - Sign of first input string in preferred form 

R7<3:0> = Sign of second input string in preferred form 


a 9 
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+ 38 ‘ «SUBTITLE ADD_PACKED = Add Two Packed Decimal Strings 
Py | 
. 435 3; Functional Description: | 
DA 461; This routine afte two packed decimal eet whose descriptors 
DA re) 3 are passed as input parameters and places their sum into another 
+ te7 3 (perhaps identical) packed decimal string. 
DA 465; At the present time, the result is placed into a 16-byte storage 
DA £96 3 area while the sum js being evaluated. This drastically reduces 
ODA 467; the number of different cases that must be dealt with as each 
p + r+ 3 3 pair of bytes in the two input strings is added. 
ODA £36 3 The signs of the two input strings have already been dealt with 
ODA 471; so this routine getreres addition in all cases, even if the original 
ODA t7¢ 3 entry was at SUBP4 or SUBP6. The cases that arrive in this routine 
ODA 473; are as follows. 
ODA 474; 
ODA 475; R2/R3 RO/R1 result 
OODA 476 ; pec mene ae ema $¢emeeseren screws t}oceereoen es 2 eece + 
QOODA 477; ' H H ; 
OODA 478; R2/R3 + RO/R1 | plus ' plus ' plus : 
QOODA 479; H : H ' 
OODA 480 3 worn arco ean ern e $¢wroanncecanrncece fp owcoenmemre caren ece Sa 
OODA 481; H : H ; 
QODA rt) 3 R2/R3 + RO/R1 } minus H minus ; minus H 
OODA 483; : H : : 
OODA 484 3 fam moeeaemene enn ane pow on nae manana we $oomewmaeneneroen se + 
OODA 485; H H ' H 
OODA 486 ; R2/R3 - RO/R1 : minus H plus H minus : 
OODA 487; ' H H ‘ 
OODA 488 3 feeee eran een eocea foe w anew wre m maw ene pew mene wera nr sana an + 
OODA 489; H : H } 
OODA 490 ; R2/R3 = RO/R1 : plus H minus H plus H 
QODA 491 ; H ' t H 
OODA 136 s ose wow wr mean ae foe mana ene mean we boawma erm we wee m ewe a Sa 
OODA 493; 
94; 
95 ; 
96 ; 
97; 
98; 
99 ; 
0; 
; 
8 : 
é 3 
8 ; 
9; 
0; 
; 
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+ i : R11 - Saved PSL (Z-bit is set, other condition codes are clear) 
DA 1 : (SP) - aoved R5, address of least significant digit of ultimate 
DA 7; destination fr ng. 
+ 18 3 4(SP) = Beginning of 20-byte buffer to hold intermediate result 
+ : Output Parameters: 
DA ¢ : The particular input operation (ADDPx or SUBPx) is completed in 
DA 3 this routine. See the routine headers for the four routines that 
ODA 4; request addition or subtraction for a List of output parameters 
DA 5 3 from this routine. 
OODA 857° 
ODA 2 $ ADD_PACKED: 
59 57 90 OODA MOVB R7,R9 ; Use sign of second string for output 
03 59 9 00DD 530 BLBC R9,10$ ; Check if sign is negative 
58 (08 «= 88 OE 331 BISB #PSL$M_N,R11 t se. SO the Saved N-bit can be set 
O0E3 5 : MARK_POINT ADD_SUB_24 
56 61 OF 8B OOE 534 10$: BICB3 #*B8000011117(R1),R6 ; Get least significant digit to R6 
Bees 535 MARK POINT ADD_S 4 
57 63 OF 8B OOE7 5 § BICB3 #*8000011117(R3),R7 ; Get least significant digit to R7 
58 Be QOEB 53 CLR ; Start the add with CARRY off 
0075 0 4454 236 BSBW ADD_PACKED_BYTE_R6_R7 ; Add the two low order digits 
OOFO 540 ; The following set of instructions computes the number of bytes in the two 
OOFO 541 ; strings and, if necessary, performs a switch so that RO and R1 always 
Bor 8 248 3; describe the shorter of the two strings. 
50 50 04 O1 EF pore 544 EXTZV #1,#4,R0,RO0 :; Convert digit count to byte count 
se 6 | US]elC SCS (itésSCéC‘irECéiF 545 EXTZV #1,#4,R2,R2 ; Do it for both strings 
52. 50 D1 OOFA 546 CMPL RO,R2 ; We want to compare the byte counts 
09 1B OOFD 547 BLEQU 20$ 3; Skip the swap if we're already correct 
56 50 7D OOFF 548 MOVG RO,R6 ; Save the longer 
~~ 2 3198 549 MOVQ R2,R0 ; Store the shorter on RO and R1 
52. 56 7D 010 550 mMOva R6,R2 3} «ee and store the longer in R2 and R3 
ie Te - Bios 22) 20S: SUBL RO,R2 ; Make R2 a difference (R2 GEQU 0) 
$108 36 3; RO now contains the number of bytes ponesney in the shorter string. 
3198 22¢ 3; R2 contains the difference in bytes between the two input strings. 
50 oS 0108 556 TSTL RO ; Does shorter string have any room? 
06 13 10P 336 BEQL 40$ 3; Skip loop if no room at all 
004D 30 O10F 25 30$: BSBW ADD_PACKED_BYTE_STRING ; Add the next two bytes together 
FA 50 FS Vg 9 SOBGTR RO,30$ ; Check for end of loop 
52 0 115 6¢ 40$: TSTL Re 3 Does Longer string have any room? 
we % at 4 o7 BEQL 70$ :; Skip next loops if all done 
0D 58 =«€9 ite 2? 50$: BLBC R8,60$ ; Life is simple if CARRY clear 
56 04 O1IC $9 CLRL R6 ; Otherwise, CARRY must propogate 
11E 68 MARK_POINT ADD_SUB_24 
57 73 QA OIE 8 MOVZBL =(R3),R7 ; So add CARRY to single string 
0041 30 0121 0 BSBW ADD_PACKED_BYTE_R6_R7 3; Use the special entry point 
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F2 5 F5 0124 71 SOBGTR Re 50$ ; Check for this string exhausted 
0 11 } 4 Me BRB 70$ ; Join common completion code 
129 74 MARK _POINT ADD_SUB_24 
7 £ 9 129 75 60$: MOVB =-(R3) ,-(R5) ; Simply move src to dst if no CARRY 
FA 52 OF ‘ ud: SOBGTR R2,60$ t se. until we're all done 
7s 3 & 0 ; 4 70$: MOVB RB,-(R5) ; Store the final CARRY 
1 80 ;+ 
61 81 ; At this point, the result has been computed. That result must be moved to 
1 ¢ ; its ultimate destination. noting whether any nonzero digits are stored 
: ? 3 so that the Z-bit will have its correct setting. 
3 282 : Input Parameters: 
0132 587: R9<7:0> = Sign of result in preferred form 
01 § 588 ; R11<3:0> = Saved condition codes 
8138 289 5 R11<31> = Indicates whether to set saved R4 to zero 
0132 «591 ; (SP) = Saved RS, high address end of destination string 
0133 893° 
0132 594 ADD_SUBTRACT_EXIT: 
55 : @ &y . Sime i ae ADDL #1,(SP),R5 ; Point RS beyond real destination 
51 18 AE 9E€ 0136 596 MOVAB 24(SP),R1 ; R1 locates the saved result 
010C 30 O13A 597 BSBW STORE _RESULT 3; Store the result and record the Z-bit 
12 5B 02 €0 He} a8 BS PSL$0_Z,R11,100$ ; Step out of Line for minus zero check 
0141 600 MARK_POINT ADD_SUB_24 
9§ 04 00 59 FO 0141 601 80S: INSV R9,#0,#4,a(5P)+ ; The sign can finally be stored 
5E 14 CO 0146 60 ADDL . ; Get rid of intermediate buffer 
03 58 1F €1 0149 60 BBC #ADD_SUB_V_ZERO_R4,R11,90$ ; Branch if 4-operand opcode 
10 AE D4 Q014D 604 CLRL 16(SP) ; Clear saved R4 to return zero 
FEAD' 31 8129 oR? 90$: BRW VAXSDECIMAL_EXIT ; Exit through common code path 
0153 607 ; If the result is negative zero, then the N-bit is cleared and the sign 
$123 one : is changed to a plus sign. 
5B 08 8A 0153 610 1008: BICB #PSL$M_N,R11 ; Clear the N-bit unconditionall 
E758 01 EO 0156 611 BBS #PSL$V_V.R11,80$ ; Do not change the sign on overflow 
59 Oc 90 Bin ol¢ MOVB #12,R9 ; Make sure that the s gn is plus 
E2 11 OSD 61 BRB 80$ 3 ee. and rejoin the exit code 


ao 
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- Add Two Bytes Containi 5-SEP-19 0:44:34 CEMULAT.SRCJVAXARITH.MAR; 1 ( 
: - SUBTITLE ADD_PACKED_BYTE = Add Two Bytes Containing Decimal Digits 
Functional Description: 


aw 
woo 
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This routine adds together two bytes containing decimal digits and 
produces a byte containing the sum that is stored in the output 
string. Each of ne input bytes is converted to a binary number 
(with a table-driven conversion), the two numbers are added, and 
the sum is converted back to two decimal digits stored in a byte. 


This routine makes no provisions for bytes that contain illegal 
decimal digits. We are using the UNPREDICTABLE statement in the 
architectural description of the decimal instructions to its fullest. 


The bytes that contain a pair of packed decimal digits can either 

exist in packed decimal strings located by R1 and R3 or they can 

be stored directly in registers. In the former case, the digits must 

be extracted from registers before they can be used in later operations 
because the sum will be used as an index register. 


| 

| 

For entry at ADD_PACKED_BYTE_STRING: | 
Input Parameters: | 
| 

| 


R1 = Address one byte beyond first Pee that is to be added 
R3 = Address one byte beyond second byte that is to be added 
oe R5 = Address one byte beyond location to store sum 


R8 - Carry from previous byte (R8 is either 0 or 1) 


R6 - Scratch 
R7 - Scratch 


Output Parameters: 
Ril - Decreased by one to point to current byte in first input string 
R3 - Decreased by one to point to current byte in second input strin 
RS - Decreased by one to point to current byte in output string 


R8 - Either 0 or 1, reflecting whether this most recent ADD resulted 
in a CARRY to the next byte. 


For entry at ADD_PACKED_BYTE_R6_R7: 
Input Parameters: 


POAOASACOOO OAD 
AOA TU 
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66 

R6 - First byte containing decimal pt pair 

R7? = Second byte containing decimal digit pair 

R5 - Address one byte beyond location to store sum 
R8 - Carry from previous byte (R8 is either 0 or 1) 


Output Parameters: 


ee ee ee ee a ee me ee ee ee ee em ee ee a ee a ee ee ee a ee a ee a a a ae ad ed 
Pw eee eel ah ah ah db db ab dh db dh a ah dh ah Ah dh ah ah a dh a dh dh dh ab dh dh a dh dh db dh ah bh dh db dh eh dh 
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645 : Implicit Input: 
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R5 = Decreased by one to point to current byte in output string 


R8 - Either 0 or 1, reflecting whether this most recent ADD resulted 
n a CARRY to the next byte. 
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Side Effects: 
R6 and R7 are modified by this routine 


RO, R2, R4, and RY (and, of course, R10 and R11) are preserved 
by this routine 


Assumptions: 

This routine makes two important assumptions. 

1. If both of the input bytes contain only legal decimal digits, then 
it is only necessary to subtract 100 at most once to put all 
possible sums in the range 0..99. That is, 

99 + 99 + 1 = 199 LSS 200 


2. The result will be checked in some way to determine whether the 
result is nonzero so that the Z=-bit can have its correct setting. 


696 ;- 
697 
698 ADD_PACKED_BYTE_STRING: 
700 MARK_POINT ADD_SUB_BSBW_24 : 

56 71 9A 701 MOVZBL <=-(R1),R6 3 Get byte from first string | 
70 MARK POINT ADD_SUB_BSBW_24 

57 73 OA 70 MOVZBL <=(R3),R7 : Get byte from second string 


705 VAXSADD_PACKED BYTE _R6_R7:: ; ASHP also uses this routine 
706 ADD_PACRED BYTE_R6_R7: 

707 movB DECIMALSPACKED_TO_BINARY_TABLECR6],- 

708 R6 3; Convert digits to binary 


PDD DPDPED SPI 


io 
oOo 


56 O000'CF46 90 


a a ee es > — 2 > an ts as —) as — 8 2 2 a 2s 2 4 ss 


felelelelelolelelolol ala ale wlelaleoleloleolelaleloleololololelolololelolole elelelol el ela elo olalo}) 


57 Q000'CF47 90 709 MOVB DECIMALSPACKED_TO_BINARY_TABLECR7J,- ; 
7 710 R7 3; Convert digits to binary 
57 56 80 7 71 ADDB R6,R7 ; Form their sum 
57 58 ~=— 80 ? 71 ADDB R8,R7 ; Add CARRY from Last step 
58 94 7 71 CLRB R : Assume no CARRY this time 
63 8F S57 91 7 714 CMPB R7,#99 ; Check for CARRY 
07 | 7 715 BLEQU 108 : Branch if within bounds | 
58 . 7 £18 MOVB #1,R8 3 propogeye CARRY to next step 
7 64 8F 82 > on SUBB ss #100, R7 : Put R7 into interval 0..99 
75 0000'CF47 6 ay 10$: MOVB DECIMALSBINARY_TO_PACKED_TABLECR7J,- 
0 : at 4 a -(R5) ; Store converted sum byte 
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- SUBTITLE SUBTRACT_PACKED = Subtract Two Packed Decimal Strings 


aon 


Dx 
Si i i i ee te ee et ee ee ee ee ee ee ee ee ee ee ee all 
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Functional Description: 


This routine takes two packed decimal strings whose descriptors 
are passed as input parameters, subtracts one string from the 
other, and places their sum into another (perhaps identical) 
packed decimal string. 


At the present time, the result is placed into a 16-byte storage 

area while the difference is being evaluated. This drastically reduces 
the number of different cases that must be dealt with as each 

pair of bytes in the two input strings is added. 


The signs of the two input strings have already been dealt with so 
this routine performs subtraction in all cases, even if the original 
entry was at ADDP4 or ADDP6. 


; Input Parameters: 


RO<4:0> - Number of digits in first input decimal string 
R1 - Address of least significant digit of first input 
decimal string (the byte containing the sign) 


: R34: 0> - Number of digits in second input decimal string 

$ R - Address of least significant digit of second input 
: decimal string (the byte containing the sign) 

: R4<h 20> - Number of digits in output decimal string 


- Address of one byte beyond least aM 1? bana digit of 
intermediate string stored on the stack 


R6<3:0> = Sign of first input string in preferred form 
R7<3:0> - Sign of second input string in preferred form 


R11 - Saved PSL (Z-bit is set, other condition codes are clear) 

(SP) - Saved R5, address of least significant digit of ultimate 
dest ination string. 

4(SP) = Beginning of 20-byte buffer to hold intermediate result 


Output Parameters: 


PUPP BS BB BEE EWI NII ronoporononnn Sma 


SOO NAMNLE WN O ODNAULS WN — O OONAULS WN -"OOONOUSwN 9A 


cov vo vVoVTVT9VT ToC C9T9T9CC0C00C0"0C00"00C0C00"0"000"00"0"0000000T090 I 


The particular input operation (ADDPx or SUBPx) is completed in 
this routine. See the routine headers for the four routines that 
request addition or subtraction for a List of output parameters 
from this routine. 


IN 
NAA O 
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; Algorithm for Choice of Sign: 


The chores of sign for the output string is not nearly so 
straightforward as it is in the case of addition. One approach that is 
often taken is to make a reasonable guess at the sign of the result. 

If the final subtraction causes a BORROW, then the choice was incorrect. 
The fign aust be changed and the result must be replaced by its tens 
complement. 


0.000000 00 09 09 09 09. 09 09 Cd CD G9 CD CD CD Od GD CD OD CD Gd CD OD CD CD GD OD Cd CD CD CD C9 CO CO CD CD CO CO CD CODCOD CD CDCD CD CDCDOCDOCDCDORODCDODOD 1 


me ee ee ed ed ed td ds = 3 8 3 8 8 = — 8 2 3 2 — = 2 ts 4 1 2 1 2 = a 
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This routine does not guess. Instead, it chooses the input string of 
the Larger absolute magnitude as the minuend for this internal 
routine and chooses its sign as the sign of the result. 

This algorithm is actually more efficient than the reasonable 

guess method and is bee ly better than a guess method that is never 
wrong. All complete bytes that are processed in the sign evaluation 
prperecererre oop are eliminated from consideration in the | 
subtraction loop, which has a higher cost per byte. 


aon 


SQV OR 


) 


The actual algorithm is as follows. (Note that both input strings have 
a} ready had Leading zeros stripped so their lengths reflect 
significant digits.) 


1. If the two strings have unequal ogy, then choose the sign of 
the string that has the longer length. 


Wr —-OODOnNOuslwnrono 97a 


OOOODMWMOMCCDOMMCON Mme 


794 
most significant byte is larger in magni tule. 


3. If the most A hb hs pony bytes test equal, then decrease the 
Lengths of each string by one byte, drop the previous most 


significant bytes, and go back to step 2. 


4. If the two strings test equal, it is not necessary to do any 
subtraction. The result is identically zero. 


Note that the key to this routine's efficiency is that high order 
+ hg that test equal in this sep are dropped from consideration in 
the more complicated subtraction loop. 


Convert digit count to byte count 


www 
moO 
ww 
mo 


; Do it for both ote ings 
R2 3; We want to compare the byte counts 
§ ; RO/R1 represent the smaller string 
; R2/R3 represent the smaller string 


04 
04 
52 


-~—<O Mmm 
> nen 


; The two input ptr ings have an equal number of bytes. Compare magnitudes to 
3; determine which string is really larger. If the two strings test equal, then 
3; skip the entire subtraction Loop. 


SUBL RO,R1,R8 
SUBL + Toma 
TST R 


STL 
BEQL 208 


0 

1 

§ 

4 

5 

: 

5 

0 SUBTACT_PACKED: 
1 EXxTZvV 
: 

4 

5 

5 

: 

0 


vw 
ooo 


Point R8 to low address end of RO/R1 
Point R9 to low address end of R2/R35 
See if both strings have zero yet 


vu 
w= 
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900 


Ww 


Still need to check low order d 
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108: CMPB (RB)+, (RO) +” : Compare most significant bytes 

; RO/R1 represent the smaller string 
BGTRU 30% 3 R2/R3 cepresent the sant tor string 
DECL Re : Keep R2 in step with R 
SOBGTR RO,10$ : 


«++ which gets decremented here 
; At this point, we have reduced both input steings to single bytes that | 
3 contain a sign ‘digit’ and may contain a digit in the high order nibble 


; if the original digit counts were nonzero. | 


2. For ocr ings of equal length, choose the sign of the string whose 
| 


no —-— 
Vine 
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MARK POINT ADD_S 

58 68 OF 88 208: B1CB3 #7B000011117(RB 
MARK _P ADD_S 
59 9 of 111,¢ 


H Get the digit from the other string 


CMPB f io : compere these ty hy 
BLSSU 408 3 aa represent the smaller string 
BGTRU 308 ; R2/R3 represent the smaller string 


; The two strings have identical magnitudes. Enter the end processing code 
3; with the intermediate result unchanged (that is, zero). 


BRW ADD _SUBTRACT_EXIT ; Join the common completion code 
The string described by RO and R1 has the Lorene magnitude. Choose its sign. 
Then swap the two string descriptors so that the main subtraction loops 
always have R2 and R3 Seer the larger string. Note that the use of 
R6 and R7 as scratch leaves R7<31:8> in an UNPREDICTABLE state. 


30$: MOVB R6,RO 


| 
R8 ; Look only at digit, ignoring sign 


FF6A 31 


; Load preferred sign into R9 


59 6 556 3 
56 28 7D RO ; Save the longer 
50 65 7D mMOVQ R2,RO ; Store the shorter on RO and R1 
i Ie. . MOVG R6,R2 3 «ee and store the honger in R2 and R3 
57 04 CLRL R7 3 Insure that R7<31:8> is zero 
03 11 BRB 50$ 3 Continue along common code path 
. 3; The string described by R2 and R3 has the larger magnitude. Choose its sign. 
40$: MOVB R7,R9 ; Load preferred sign into R9 
52 C2 50$: SUBL RO,R2 ; Make R2 a difference (R2 GEQU 0) 
03 59 + 6 BLBC R9,60$ ; Check if sign is negative 
58 8 BISB #PSLSM_N,R11 3 ee. SO the saved N-bit can be set 


MMMMOOVTOCVCVVC TCV IGTAAAAAOOOAOOOOOOOO@OWWOew «~«— 


WW OUUIW SMO COCO & & HM DW W COC COCOS: & — FN W COCO COCO COCO OOF DB BAAD 


MARK_POINT ADD_SUB_24 
60$: BICB3 #*B000011117(R1),R6 : Get least significant digit to R6 
MARK_POINT 


ADD 
BICB3 #*8000011117(R3),R7 : Get least significant digit to R7 
CL ; Start subtracting with BORROW off 
BSBW SUB_PACKED_BYTE_R6_R7 ; Subtract the two low order digits 


; RO contains the number of bytes remaining in the smaller string 
3; R2 contains the difference in bytes between the two input strings 


7 
TSTL RO ; Does smaller string have any room? 

1 BEQL 80$ ; Skip loop if no room at all 

: 70$: BSB SUB_PACKED_BYTE_STRING ; Subtract the next two bytes 

4 RO,70$ L 

§ 80$ TSTL 
3 

2 


WW 
SOBGTR 7 : Check for end of loop 
Does one of the strings have more? 


56 61 OF 8B 
57 63 


R2 
BEQL 110$ Skip next loops if all done 
90$: BLBC R8,100$ ; Life is simple if BORROW clear 


EDI a a a kn ak ht a 8 8 oe 
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FUN TS 


CLRL R6 ; Otherwise, BORROW must propogate 
MARK_POINT ADD_SUB_24 


SS ——— aed 
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Subtract Two Packed De 5-SEP-1984 00:44: EMULAT.SRCJVAXARITH.MAR; 1 (11) 


MOVZBL a. ; So subtract BORROW g-  gingte string 


-. FT A 
01 30 BSBW = SUB ere _BYTE_R6_R7  ; Use the aeeetol entry 
a 3 F SOBGTR R2,90$ : Check for this string Serecsted 
11 BRB 1 65" : Join common completion code 
MARK _POINT ADD_SUB_24 
100$: MOVB™ =(R3),-(R5)~ ; Simply move src to dst if no BORROW 
SOBGTR R2, 106$ 3 ee. until we're all done 


110$: 


seg tteeeeeeee BEGIN TEMP teenenenee 


ti: THE FOLLOWING HALT INSTRUCTION SHOULD BE REPLACED WITH THE CORRECT 
$i: ABORT CODE. 

ti: THE HALT IS SIMILAR TO THE 

ih: MICROCODE CANNOT GET HERE 


333 ERRORS THAT OTHER IMPLEMENTATIONS USE. 


tstl r8 : If BORROW is set here, we blew it 
beql 1208 : ecenen out if 0 
halt 3; This will cause an OPCDEC exception 


oul 
—co 
ee 


dS 
13 
00 
120$: 


tee 
: weeeerererte END TEMP teeeaannane 


BRW ADD_SUBTRACT_EXIT ; Join common completion code 
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v04 


1984 01:33:44 Page 21 
- Subtract Two Bytes Con met 7 Bhi 23:32 YEMULAT. SREIVAXARITH.MARS 1 . (ad) 


- SUBTITLE SUB_PACKED_BYTE = Subtract Two Bytes Containing Decimal Digi 
Functional Description: 


This routine takes as input two bytes containing decimal digits and 
produces a b ? containing their difference. This -esult is stored in 
the output string. Each of the ingus oytes is gonyerves to a binary 

th on), the first number is 
subtracted from the second, and the difference is converted back to 
two decimal digits stored in a byte. 


number (w a table-driven convers 


This routine makes no provisions for bytes that contain illegal 
decimal digits. We are us ng the UNPREDICTABLE statement in the 
architectural description of the decimal instructions to its fullest. 


The bytes that contain a pair of packed decimal digits can either 

exist in packed decimal strings located by R1 and R3 or they can 

be stored directly in registers. In the former case, the digits must 

be extracted from registers before they can be used in later operations 


1 

K 

10 

1D H 

10 : 

1D ; 

10 3 

1D 3 

1D 3 

1D : 

10 3 

10 3 

1D 3 

1D ; 

1D 5 

1D 3 

1D 3 

10 3 

1D : 

B 3: 

1 Be ; because the difference will be used as an index register. 

1p 94s : For entry at SUB_PACKED_BYTE_STRING: 

1 $47 : Input Parameters: 

1D 949 : Ri = Address one byte beyond byte containing subtrahend 

1D 950 ; R3 = Address one byte beyond byte containing minuend 

4 33) 5 R5 = Address one byte beyond location to store difference 

4 $58 : R8 = BORROW from previous byte (R8 is either 0 or 1) 

1p 935 : Implicit Input: 

1D «957: R6 = Scratch 

19 228 3 R7 = Scratch 

1 44 : Output Parameters: 

1D 96 : R1 = Decreased by one to point to current byte 

1D 963 ; n subtrahend string 

1D 964 ; R3 - Decreased by one to point to current byte 

1D 965 ; n minuend strin 

1D 208 3 R5 - Decreased by one to point to current byte 

19 $ $ in difference string 

1D 569 ; R8 - Either 0 or 1, reflecting whether this most recent 

19 af § subtraction resulted in a BORROW from the next byte. 

1p 97¢ : For entry at SUB_PACKED_BYTE_R6_R7: 

1 afe : Input Parameters: | 

1D 378 : R6<7:0> = Byte containing decimal digit pair for subtrahend 

1b 977 ; R6<3i:8> - nbz 

1D 4 3 R?7<7:0> - axe containing decimal digit pair for minuend 

1D 979; R7<31:8> - MBZ 

1D 980; 
| 


<< 
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VOee 000 SUB_PACKED_BYTE = Subtract Two Bytes Con 5-SEP-1984 bb ita: $6 EMULAT.SRCJVAXARITH.MAR; 1 ’ (3) 
} 1 : R5 - Address one byte beyond location to store difference 
1 : ; R8 = BORROW from euptrece ton of previous byte 
: ¢ : (R8 is either 0 or 1) 
: : § ; Output Parameters: 
1D «988: RS = Decreased by one to point to current byte 
: 9 5 3 in difference serine” r 
1 391 F R8 - Either 0 or 1, reflecting whether this most recent 
! $38 : subtraction resulted in a BORROW from the next byte. 
934 : Side Effects: 
936 : R6 and R7 are modified by this routine 
1 998 : RO, R2, R4, and RI (and, of course, R10 and R11) are preserved 
: +4 3 by this routine 
001 ; Assumptions: 
; 098 : This routine makes two important assumptions. 
1D 1005 ; 1. If both of the input bytes contain only legal decimal digits, then 
1 4 3 it is only necessary to add 100 at most once to put all 
: 44 3 possible differences in the range 0..99. That is, 
09 ; 0-99 -1 = -100 
1 11: The result will be checked in some way to determine whether the 
; ; $ result is nonzero so that the Z-bit can have its correct setting. 
1 14° 
: 0 SUB_PACKED_BYTE_STRING: 
1 01 MARK_POINT ADD_SUB_BSBW_24 
56 71 OA 1 1 MOVZBL =(R1),R6 3 Get byte from first string 
1 MARK _POINT ADD_SUB_BSBW_24 
57 73 9A MOVZBL -(R3),R7 ; Get byte from second string 


SUB_PACKED _BYTE_R6_R7: 
move DECIMALSPACKED_TO_BINARY_TABLECR6),- 
3_Conv rt digits to binary 
MOVB DECIMALSPACKED_TO_BINARY_TABLECR7]),- 
R7 ; Convert digits to binary 
SUBB R6,R7 : Form their difference 


56 O000°CF46 90 
57 Q000°CF47 90 


DNV NOOO OSC CUT VT TV VTVCTVVVVVOVCCUVCVTVCUVCVUCOCCUCCCOCVCVCOVCOVCUVCVUVUCUTUTCSO 


57 

57 $ 83 SUBB R8,R7 Include BORROW from Last step 
Qe | BLSS los : Branch if need to BORROW 
: 94 CLRB Rp No BORROW next time 
0 11 BRB 20$ Join common exit code 


R7 ; Put R7 into interval 0..99 
§ ; Propogate BORROW to next step 


20$: MOVB DECIMALSBINARY_TO_PACKED_TABLECR7),- 
-(R5) 3; Store converted sum byte 


Tse O90 


75 O000'°CF47 90 


baal 


10$: ADDB #100 
MOVB #i1,R 


SOOOSOOOSOSOSSOOOCOOOOOCOOSOOSOSOSOOOSCOOOOCOSCOOSOOCSOOCOOCOOOOOOOOOOOOOOOO 


ern 
Coron 
a a ee ee ee ee ee ee eee ee ee ee 


a YF 
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Ye aaa SUSAPACKED BYTE "Subtract fwo Bytes ton ‘S-SEP-108¢ Qd:4e:4¢ ERMOLAT.SRESWAXARION.maR;1 29° (93) 
05 0248 1038 RSB 
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vee ts STORE RESULT = Store Decimal String ~$Eb= 1382 8b i224 EMULAT.SRCJVAXARITH.MAR; 1 at F | 
r } "4 ie - SUBTITLE STORE_RESULT = Sture Decimal String 
4 4 ; Functional Description: | 
49 1 44 : This routine takes packed dec inal string that typically contains 
49 10465 ; the result qf n arithmetic operation and stores it in another 
$3 1 rf 3 decimal string whose eooes ptor is specified as an input parameter | 
rh : : : to the original arithmetic operation. 
49 1 rk : The string is stored from the high address end (least significant 
49 1050 ; digits) to the low address end (most significant digits). This order 
rh ! ; allows all of the special cases to be handled in the simplest fashion. 
49 33 ; Input Parameters: 
49 1055 : R1 - Address one byte beyond oy address end gf input string 
re : 2$ ; (Note that this string must be at least 17 bytes long.) 
49 1988 ; R4<4:0> - Number of digits in ultimate destination 
re ' re : R5 - Address one byte beyond destination string 
rh 1961 : R11 - Contains saved condition codes 
0 rh 1068 : Implicit Input: 
6 49 1985 : The input string must be at least 17 bytes long to contain a potential 
0249 1 96 3 carry out of the highest digit when do ng an add of two large numbers. 
8 49 1067 ; This carry out of the last byte will be detected and reported as a 
re 1908 3 decimal overflow, either as an exception or simply by setting the V-bit. 
49 19%6 : The least significant digit (highest addressed byte) cannot contain a 
0 rh 8 : sign digit because that would cause the Z-bit to be incorrectly cleared. 
0 rh 1078 : Output Parameters: 
0249 1075 : R11<PSL$V_Z> - Cleared if a nonzero digit is stored in output string 
0249 1376 3 R11<PSL$V_V> = Set if a nonzero digit is detected after the output 
0 rt 4 : 4 ; string is exhausted 
49 1 8 : A portion of the result (dictated by the size of R4 on input) is 
rh 4 } , 3 moved to the destination string. 
49 1082 © 
rh 1 § STORE _RESULT: 
54 49 1084 INCL RG 3; Want number of ‘‘complete’’ bytes in | 
50 54 FF 8F 48 1085 ASHL #-1,R4,R0 3 output string 
0B 1 ! § BEQL 30$ 3; Skip first loop if none 
1 3 MARK_POINT ADD. SUB_BSBW_24 
, aa. oe 1 10$: MOVB ~(R1),-(R5) 3; Move the next complete byte 
ip 55 1090 BEQL 3; Check whether to clear Z-bit 
58 4 BA 1091 BICB #PSL$M_7Z,R11 ; Clear Z-bit if nonzero 
F550 =sF 5 5 | 35 20$: SOBGTR RO,10$ 3; Keep going? 
10 54 =€9 1094 30$: BLBC R4,50$ ; Was original R4 odd? Branch if yes | 
69 1095 MARK_POINT ADD_SUB_BSBW gf 
75 71 FO BF 8B 0260 1096 BICB3 #*B111100007-(RT) ,-(R5) ; If R4 was even, store half a byte | 
j 
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“SEP-1984 01:33:44 VAX/VMS M v04-00 
STORE_RESULT = Store Decimal String g SE 3b §bi2a:44 f REQ 


P= EMULAT.SRCJVAXARITH.MAR; 1 


oan ; 1097 BEQL $ 3; Need to check for zero here, too 
58 4 BA 1 38 BICB #PSLSM_7,R11 3; Clear Z-bit if nonzero 
A 109 MARK_POINT ADD_SUB_BSBW_24 
61 FO 8F 3 A 1 40$: BITB €oBt11 10000, Cn) : If high order nibble is nonzero, 
1 1 8 i 1 BNEQ 0 3 eee then overflow has occurred 
70 «#11 7 ; The entire destination has been stored. We must now check whether any of 
70 1106 ; the remainin input string is nonzero and set the V-bit if nonzero is 
70 1105 ; detected. Note that at least one oree of the output string has been examined 
f i $ : in all cases already. This makes the next byte count calculation correct. 
4 07 70 1108 50$: DEC RG 3; Restore R4 to its original self 
50 4 046 t ff 7 1199 PxTe #1,44,R4,R0 ; Extract a byte count . 
0 10 0 3 A 1119 SUBB RO.#16,RO ; Loop count is 16 minus byte count 
78 1116 ; Note that the loop count can never be zero because we are testing a 17-byte 
A a 3; string and the largest output string can be 16 bytes long. 
78 1115 MARK_POINT ADD_SUB_BSBW_24 
71 395 78 1138 60$: TSTB =(R1) 3; Check next byte for nonzero 
9 ig 70 111 BNEQ 70$ : Nonzero means overflow has occurred 
F9 50. CéF 8 ze 1133 SOBGTR RO,60$ 3; Check for end of this loop 
05 8 e3 i RSB ; This is return path for no overflow 
5B 02 88 0283 11 ¢ 70$: BISB #PSLSM_V,R11 ; Indicate that overflow has occurred 
05 0286 11 RSB 3 ee. and return to the caller 
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r 
EMULAT.SRCJVAXARITH.MAR; 1 
- SUBTITLE VAXSMULP = Multiply Packed 


Functional Description: 


The multiplicand string specified by the mul tipticend length and 
multiplicand address operands is multiplied by the multiplier string 
specified by the multiplier length and multiplier address operands. The 
product string specified by the product length and product address 
operands is replaced by the result. 


Parameters: 

RO = mulrlen.rw Number of digits in multiplier string 
R1 = mulraddr.ab Address of multiplier ete ing 

R2 - muldlen.rw Number of digits in multiplicand string 
RS = muldaddr.ab Address of multiplicand strin 

R4 = prodlen.rw Number of digits in product string 

RS = prodaddr.ab Address of product string 


Output Parameters: 


r] 
oO 
i] 


= 0 
Address of the byte containing the most significant digit of 
~ multiplier string 


»@ 
= 
" 


Address of the byte containing the most significant digit of 
is multiplicand string 


»@ 
ea 
“uu 


Address of the byte gale the most significant digit of 
the string containing the product 


Condition Codes: 


N <= product string LSS 0 
2 <= product etrioy EQL 0 
V <= decimal overflow 


C <- 0 


Register Usage: 


PIPPIN IPIDIMNIMNIDININININIPNIPYNIPYNIPINININININININININININY[PINIPUNIPIPYAIPUAIPIAIPPIPIPIPIPUNPIPOPINIPUNPYD 
0.09.09 09 Gd Gd Cd Gd Gd Cd Cd Cd GD Gd Gd Cd Cd OD OD CD Cd OD Od 0D OD OD CD Od 0D GD CD OD OD OD CD OD OD CS CD CD CD CD CD CD OP CD CD C9 G9 OD CO CD CD OD CD CDCD 
i i i Bi Be i Bt te DD te fe tt tt tt te i ee te ee i ee ee 
a ee ee a ee ee et me ee ee ee a ee ae a a 
DOOD INN NSD DDD DEDEDE DDS TUT Be Be ES BB BB BWI I nnrnnuny) 
=O ODNOAO UNE WIN SO OD NOME WIN OS OD NA UE WIN 0 OD NAME WO OONOA UE WIN OOONOu 


COOOCOCOOOOCOCOCOCOOCOOCOCOCOOOOOOOOCOCOCOOCOOOOOOOOOOoOOo 


This routine uses all of the general registers. The condition codes 
are computed at the end of the instruction as the final result is 
eteres n the product string. R11 is used to record the condition 
codes. 


This routine uses a large amount of stack agece to allow etereay of 
intermediate results in a convenient form. Specifically, each digit 
pair of the Longer input string is stored in binary in a lLongword on 
the stack. In addition, 32 longwords are set aside to hold the product 
sntergedsate result. Each longword contains a binary number between 
an 


After the multiplication is complete. Each longword is removed from 
the stack, converted to a packed decimal pair, and stored in the 
output string. Any nonzero cells remaining on the stack after the 


0 v04-00 p 6 | 
age 98, | 
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v04- VAXSMULP = Multiply Packed -SEP-1984 00:44:34 CEMULAT.SRCJVAXARITH.MAR;1 (14) 
Stee ¢ ; output string has been completely filled are the indication of decimal 
f 1 ? : overflow. 
7 1185: The purpose of this method of storage is to avoid decimal/binary or 
y 6(«Y § : even byte/longword conversions during the calculation of intermediate 
f 1 : results. 
7 1189 ; 2. Trailing zeros are removed from the Larger string. ALL zeros in 
7 1190 ; the shorter string are eliminated in the sense that no arithmetic 
0287 1191; is performed. The output array pointer is simply advanced to point 
0287 1136 3 to the next higher array element. 
8 sf 1193 ;:- 
7 1194 
0287 1195 VAXSMULP:: 
OFFF 8F BB 3 rH 1138 PUSHR #*°M<RO,R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,R11> ; Save the lot 
0388 1198 ESTABLISH HANDLER - ; Store address of access 
th) 113 ARITH_ACCVIO 3 violation handler 
0390 1201 ROPRAND_CHECK R4 ; Insure that R4 is LEQU 31 | 
0298 1208 ROPRAND_CHECK : Insure that R2 is LEQU 31 
bsAs 1204 MARK_POINT MULP_BSBW_0 / 
FDSA' 30 BSA? 1599 BSBW DECIMALSSTRIP_ZEROS_R2_R3 ; Strip high order zeros from R2/R3 
02A6 1207 ROPRAND_CHECK : Insure that RO is LEQU 31 
O2AE 1208 MARK_POINT MULP_BSBW_0 
FD4F* 30 BgAt 159) DECIMALSSTRIP_ZEROS_RO_R1 ; Strip high order zeros from RO/R1 
ss U©R. SR TH OU 6561 1$1} EXTZV #1,#4,R0,RO0 ; Convert digit count to byte count 
50 06 bSe6 \si¢ INCL RO : Include least significant digit 
52 52 04 O1 EF 0288 1214 EXTZV -#1,#4,R2,R2 ; Convert digit count to byte count 
52.6 pgB0 1$1? INCL R2 : Include least significant digit 
S250 = «(1 $58 1217 CMPL RO,R2 : See which string is larger ; 
08 iA O2C2 1218 BGTRU 3$ ; R2/R3 describes the penper string 
58 52 7D O2C4 1219 MOVa R2,R8 ; R8 and RY describe the longer oer ing 
7E 50 7D Q2C7 \§ 0 mMOVQ RO,-(SP) ; Shorter string descriptor also save 
06 «11 Beta 1s 1 BRB 6$ 
= 2 PP Bsce 1 : 3$: MOVQ RO,R8 ; R8 and RO describe the longer strin 
oe ee cf ’ 3 mova R2,-(SP) ; Shorter string descriptor also save 
020 ' § ; Create space for the output array on the stack (32 longwords of zeros) 
50 08 00 : \ 3 6$: MOVL #8,RO ; Eight pairs of quadwords | 
ae D5 1250 10$: CLRQ -(SP) 3; Clear one pair 
mm: D7? 1231 CLRQ -(SP) 3 ee. and another 
F9 50—osF 5 4 : ¢ SOBGTR RO,10$ ; Do all eight pairs | 
57 SE 00 be } ; MOVL SP,R7 ; Store beginning of output array in R7 | 
DF 1 § 3; The banger input array will be stored on the stack as an array of 
DF 1 ; longwords., Each array element contains a number between 0 and 99 
DF 1258 ; representing a pair of digits in the original packed decimal string. 


VAXSDECIMAL ARITHMETIC = VAX- 
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D 
iply Packed “SEP-19 0:44:34 CEMULAT.SRCJVAXARITH.MAR; 1 
Because the units digit ts stored with the sign in packed decimal format, 


it is necessary to s the number as we store it. This is accomplished by 
multiplying the number by ten. 


The longer array is described by R8 (byte count) and R9 (address of most 
significant digit pair). 


a> Soe, oes ADOL3 = R9,R8B,R5 ; Point RS beyond sign digit 
54 53 00 MOVL RB,R4 : R4 contains the loop count 
3 An array of longwords is allocated on the stack. R3 starts out pointing 
; at the longword beyond the top of the stack. The first remainder, guaranteed 
; to be zero, is ‘'stored’’ here. The rest of the digit pairs are stored safely 
; below the top of the stack. 
5 5 CE MNEGL R8,R3 ; Stack grows toward lower addresses 
E 6643 »D MOVA (SP)ER3}.,SP : Allocate the space 
33...§ 06 =«CC SUBL #4,SP,R : Point R3 at next lower longword 


MARK_POINT MULP_R8 
20$: MOVZBL -(RS),R1 et next digit pair 


ee 3: G6 
51 O000'CF41 9A MOVZBL DECIMAL SPACKED_T0_BINARY TABLECR1J,- 


3 Convert digits to binary 
50 52 51 OA EMUL #10,R1,R2,R0 3 ey eipty by 10 
83 52 50 00000064 8F 78 EDIV. #106,R6,R2,(R3)+ : Divide by 100 
£6 54 OFS 08 SOBGTR R4,20$ 
MOVL R2,(R3) ; Store final quotient 
MOVL ’ ; Remember array address in R9 
PUSHAL (SP)CR8) 3; Store start of fixed size area 


; Check for trailing zeros in the input array stored on the stack. If any are 
3; present, they are removed and the product array is adjusted accordingly. 


30$: TSTL (R9)+ 
BNEQ 40$ 


ADDL #4,R7 
SOBGTR R8,30$ 


:; Is next number zero? 

: Leave loop if nonzero 

3; Advance output pointer to next element 
; Keep going 


; If we drop through the loop, then the entire input array is zero. There is 
3; no need to perform any arithmetic because the product will be zero (and the 
3; output array on the stack starts out as zero). The only remaining work is 

; to store the result in the output string and set the condition codes. 


=O OONAUE WN OOD NAME WN SO ODNAUE WN OOCONAUEWN0O0 7a 


BRB 70$ ; Exit to end processing 


Now multiply the input array by each successive digit pair. In order to 
allow R10 to continue to locate ARITH_ACCVIO while we execute this Loop, 
is necessary to perform a small amount of register juggling. In essence, 
R8 and RY switch the identity of the string that they describe. 


40$: SUBL #4,R9 peadivat i.put array pointer 
MOVQ R8,-(SP) Save R3/R9 descriptor_on stack 
MOVL 8 (SP) -R8 Point R8 at start of 32-longword array 
MOVa <3224>(R8B) RB 


Get Coser iptor that follows that array 


ADDL2 R8,R9 Point R9 beyond sign byte 


WOO OO O0090900 09 09 0909 G9 00 SI NIN NNN SDS DDD AAO I EEE EEE OO 


MEW OOOO 
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vive tit VAXSMULP = Multiply Packed -SEP-1984 ieee EMULAT.SRCJVAXARITH.MAR; 1 . d 
53 87 ODE 1 36 50$: MOVAL (R7)+,R3 ; Output array address to R3 
129 MARK _POINT MULP_AT_SP 
51 79 9A 1 4 MOVZBL <=(R9),R1 :_Next digit pair to R1 
56 O000'CF41 9A 109 MOVZBL DECIMALSPACKED_TO_BINARY.TABLELR11,- 
te ? R6 ; Convert digits to binary 
06 13 1 BEQL 60$ ; Skip the work if zero 
54 fe. 40 1 ¢ MOVQ (SP) ,R4 ; Input array descriptor to R4/R5 
Qi 4 3 43 1 BSBW EXTEND STRING_MULTIPLY ; Do the wor 
e958 SOF 8 : : 60$: SOBGTR R8,50$ ; Any more multiplier digits? 
5E 08 CO rts ! 5 ADDL #8,SP ; Discard saved long string descriptor 
5E 6E€ 00 rt: } 38 70$: MOVL (SP) ,SP ; Remove input array from stack 
4F 1310 ; At this point, the produce string is located in a 32-longword array on 
34F 1311 ; the top of the stack. Each longword corresponds to a pair of digits in 
O34F 1 \¢ ; the output ering. As digits are removed from the stack, they are checked 
0O34F 1313 ; for nonzero to obtain the correct setting of the Z-bit. After the output 
O34F 1314; string has been filled, the remainder of the product string is removed from 
tat 1312 ; the stack. If a nonzero result is detected at this stage, the V-bit is set. 
59 20 D0 O34F 1319 MOVL #32,R9 ; Set up array counter 
54 0098 CE 7D 039¢ 1318 MOVO << <$2e4> + - ; Skip over 32-longword array 
0357 1319 <2t4> + = ; and saved string descriptor 
0357 1320 <4*4> >(SP),R4 ; to retrieve original R4 and R5 


] 
' 


9. 
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» SUBTITLE Common Exit Path for VAXSMULP and VAX$DIVP 
The code for VAXSMULP and VAXSDIVP merges at this point. The result is stored | 
in an array of longwords at the cop of the stack. The size of this array is 
stored in R9. The original R4& and R5 heve been retrieved from the stack. 


Input Parameters: 
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R4 = Contains byte count of destination greene in R4 <1:4> 
R5 = Address of most significant digit of destination string 
R9 = Count of longwords in result array on stack 
Contents of result array 
Implicit Input: 


Signs of two input factors (multiplier and multiplicand or 
divisor and dividend) 


z Se Se Ge Ge Ge Ge Ge Ge Se Ge Ge Ge Ge Ge Ge Ge Se Se 
4 * 


ULTIPLY_DIVIDE_EXIT: 
OVPSL “R1 


A. AANA AAI AIA AIAINI AIAN WI NININIAI AAAI AAAI AIAN AAAI AAI AANA 


WN ODNAUE WIN 9 OD NA UE WIN O OD NAME WN OOD NAME WO ODN UNE Un 


MARK_POINT MULP_DIVP_R9 
110$: #MOVB DECIMALSBINARY_TO_PACKED_TABLECR1),- 
-(R5) ; Store converted sum byte 
pec. R9 : One less element on the stack 


1168 
BRB 140$ 


75 0000'CF41 90 


DPR VAS LLLP LILIAN ANNAN ES BES ES EEE EE AANA AA AAA IOI II 
uw 


Exit loop if result array exhausted 
; Perform overflow check 


ee me ed a a ad bs — bb — — 2 bb bs sb a bs 2b td — 2s — 2) — _ _ _ 


5 
5 
‘ 
5 
035 
8 5 
5 
ne 
035 
8 5 
5B OC 5 A 3; Get current PSL 
5B 046 00 04 FO O35 INSV #PSL$M_7,40,44,R11 3; Clear ail codes except Z-bit 
035 ESTABLISH_HANDCER - ; Store address of access 
8 5 ABITHLAGCYIO ; violation handler again 
ys 2: Ca... SS 6 EXTZV #1,#4,R4,R :; Excess byte count to R3 
13 036 BEaL «=s«125$ : Skip to single digit code 
2 tee C1 036 ADOL3) R3,R5,R7 ; Remember address of sign byte 
55 57 O11 C71 8 ° 3 ADDL3 = #1,R7,R5 3; Point RS beyond end of product string 
51 3 0 037 3 80$: MOVL (SP)+,R1 : Remove next value from stack | 
0 13 037 BEQL : Do not clear Z-bit if zero 
58 (04 ~=SOBA 0 , : BICB2 #PSL$M_Z,R11 : Clear Z-bit | 
037 ; MARK_POINT MULP_DIVP_R9 
75 0000°CF41 90 037 90$: MOVB DECIMALSBINARY_TO_PACKED_TABLECR1],- 
038 3 -(R5) ; Store converted sum byte 
59 D7 8 8 DECL R 3 One less element on the stack 
ae H BLEQ 116$ 3; Exit loop if result array exhausted 
€B 53 =F 5 ; SOBGTR R3,80$ ; Keep going? 
2254 «2&9 8 100$:  BLBC R4,120$ : Different for even digit count | 
0 ; The output string consists of an odd number of digits. A complete eigte 
3 pair can be stored in the most significant (lowest addressed) byte o | 
3; the product string. 
51 a D MOVL (SP)+,R1 ; Remove next value from stack 
0 1 BEQL 1 : Do not clear Z-bit if zero 
5B 04 BA BICB2 #PSL$M_Z,R11 ; Clear Z-bit 
78 


59 «(7 
7 15 
8 11 
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9E 1379 
of 1380 ; This a executes if the result array has fewer elements than the output 
9 1 ; string. The remaining bytes in the output string are filled with zeros. 
3 ! : There is no need for an overflow check. 
9E 1384 MARK _POINT MULP_DIVP_8 
B 96 9E 1385 114$: R =(R 3 Store another zero byte 
FB 5 Fu - } $ 116$:  SOBGEQ R3,114$ 3; Any more room in output string 
38 . i g BRB 150$ ; Determine sign of result 
8 AS 1 39 :; This code path is used in the case where the output digit count is 0 or 1. 
. : ; ; RS must be advanced 
57 22 dO 6 AS 1 38 125$:  MOVL R5,R7 ; Remember gesrese of output sign byte 
5 6006 AB 1394 INCL R5 3; Advance R5 so common code can be used 
DB Séd'YI 8 m : 2 BRB 100$ ; Join common code path 
beac 1 39 ; The output string consists of an even number of ongits. Only the low order 
O3AC 1398 ; nibble is stored in the most significant (lowest addresses) byte. A ero is 
O3AC 1399 ; stored in the niga order nibble. If the high order digit would have been 
beac 1400 ; nonzero, the V-bit is set and the overflow check is bypassed because there 
3AC 1401 ; are faster ways to clean the stack if we do not have to check for nonzero 
8 a 1258 3 at the same time. 
51 8€ DO OSAC 1404 120$: MOVL (SP)+,R1 3; Remove next value from stack 
51 QO00'CF41 90 O3AF 1405 MOVB DECIMALSBINARY_TO_PACKED_TABLECR1j,- 
tt 1406 ; Obtain converted sum byte 
385 1407 MARK_POINT MULP_DIVP_R9 
75 51 FO 8F 58 p38 1408 BICBS #*xXF0,R1,-(R5) i Store byte, clearing high order nidble 
_— or 1409 EQL 3; Do not clear Z-bit if zero 
5B 04 8A OQ38C 1410 BICB2 #PSL$M_7,R11 : Clear Z-bit 
51 FO 8F 93 Q3BF 1411 130$: BITB #*xFO,R1 :; Is high order nibble nonzero? 
96 Ig 03¢3 1716 BNEQ 3; Yes, go set overflow bit 
9 od? 0 ce 141 DECL R9 : One less element on the stack 
D7 «15 O3C7 1414 BLEQ 116$ 3; Exit loop if result array exhausted 
0B 11 2 34 131? BRB 140$ ; Check rest of result array for nonzero 
03¢B 1419 3; If we detect overflow, we need to adjust R9 to reflect the nonzero longword 
03CB 1418 ; removed from the stack before we enter the next code block that sets the 
. 101? 3; V-bit and cleans off the stack based on the contents of R9. 
59 07 8 1? 1 133$: DECL R9 3; One more longword removed from stack 
CD 14 : 3; A nonzero digit has been discovered in a position that cannot be stored in 
CD 1424 ; the output string, Set the V-bit, remove the rest of the product array from 
CD 1425 ; the stack, and join the exit processing in the code that determines the sign 
2 1 § 3; of the product. 
58 0¢ 88 CD 14 3 135$: BISB #PSL$M_V,R11 : Set the overflow bit 
5E 664 DE DO 14 MOVAL SP RCRI Ms : Clean off repeining product string 
07 #1 1? ? RB 1508 3; Go to code that determines the sign 
44 14 § 3; The remainder of the product array must be removed from the stack. A nonzero 
D6 1435 ; result causes the V-bit to be set and the rest of the loop to be skipped. 
4 1? : ; Note that there is always a nonzero loop count remaining at this point. 
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of ) D6 14 § 140$: TSTL (SP)+ ; Is next longword zero? 
F 1 D8 14 N 153$ ; No, leave loop 
F9 59—sdéF + 1? § SOBGTR R9,140$ | 
DD 1440 ; The final product string has been stored and the V- and Z-bits have their 
DD 1441 ; correct settings. The s gn of the product must be determined from the 
DD 1208 3; signs of the two input strings. Sypes te signs produce a negative product. 
4 1227 ; Same signs (in any representation) produce a plus sign in the output string. | 
E 8 C8 CD 1445 150$: ADOL #8,SP ; Discard saved string descriptor 
6 C 4 E 1468 MOVL #12,R6 ; Assume final result is positive 
0 6€ 7D O3E3 144 mova (SP RO : Retrieve original RO/R1 pair 
50 50 046 OF EF O3E 1448 EXTZV #1,#4,R0,RO Get byte count for first input string 
1 0 co EB 144 L ‘ ; Point R1 to byte containing sign 
EE 1450 MARK_POINT BuLP DIVP 8 
50 61 FO 8F 88 e Vee] BICB3 #*B811110000,TR1).R : RO contains the sign ‘digit"’ 
3 1288 CASE 88 TYPE=B,LIMIT=#10,<- ; Dispatch on sign digit 
FS 1454 Os,- : 10 => sign is "+" 
F3 1455 10$,- ; 11 => sign is ‘'="' 
O3F 136 20$,- 3 \¢ => sign is ‘+ 
O3F3S 145 10$,- 3 13 => sign is ‘'="' 
O3F3 1458 20$,- ; 14 => sign is ‘'+"' 
8 F3 1459 20$,- : 15 => sign is ‘+’ 
O403 1461 ; 
54 01 DO 040 1066 210$:  MOVL #1,R4 ; Count a minus sign | 
02 11 0208 1007 BRB 236$ 3; Now check second input sign 
54 04 408 1692 220$: CLRL R4 ; No real minus signs so far 
52 08 AE 7D Q60A 1469 230$: mova 8(SP),R ; Retrieve original R2/R3 pair 
52 52 06 ran rh 1468 EXTZV #1,#4,R2,R2 ; Get byte count for second input string 
53 2 CO 041 1297 ADDL e : Point R3 to byte containing sign 
D416 1470 MARK_POINT MULP_DIVP_0 Se 
52 63 FO 8F 88 3 | 4 BICB3 #*B811110000,TR3) ,R2 : R2 contains the sign ‘‘digit 
41B 1298 CASE Re TYPE=B,LIMIT=#10,<- ; Dispatch on sign digit 
418 1474 6s,- : 10 => sign is +" 
041B 1475 40$,- : 11 => sign is "'-" 
418 1676 50$,- 3 \¢ => sign is ‘'+"' 
418 147 40$,- ; 13 => sign is ‘'-"' 
418 1028 50$,- ; 14 => sign is ‘'+"" 
418 147 50$,- : 15 => sign is ‘+ 
ab ue ; 
42B 14 § soos: INCL R4 ; Remember that sign was minus 
42D 14 50$:  BLBC R4,260$ ; Even parity indicates positive result 
430 1484 BBS #PSL$V_2,R11,270$ ; Step out of Line for minus zero check 
434 1485 BISB #PSLSM_N,R11 ; Set N-bit in saved PSW 
" i 1? § 255$: INCL 6 3 Change sign to minus 
439 14 g MARK_POINT MULP_DIVP_0 ; 
67 9 1489 260$:  INSV R6,#0,44, (R77 : Store sign in result string 
43E 1490 CLRL 3-16 ¢ SPS : Set saved R4 to zero 
vey 1433 BRW VAXSDECIMAL_EXIT ; Join common exit code 


H 11 
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493 ; If the result is negative zero, then it must be changed te gos itive zero 
133 3; unless overt a occurred, in which case, the s gn is left as negative 
495 ; but the N-bit is clear. 


1496 
wm 8 1497 2708: BBS —- #PSLSV_V,R11,2558 ; Make sign negative if overfl 
ei 1498 BRB it eae t Sign will be positive 
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Functional Description: 
This routine multiplies an array of numbers (each array element LEQU 


by a number (also LEQU 99). The resulting product gcray is added 
to another array, each of whose elements is also LEQU 99. 


R3 - Pointer to output array 
R4 = Input array size 
RS - input array address 
R6 = Multiplier 
Output Parameters: 
None 


Implicit Output: 


| 
| 
| 
| 
Input Parameters: | 
| 


The output array is altered. 

An intermediate product array is a by multiplying each input 

array element by the multiplier. Each product array element is then 
added to the corresponding output array element. 
Side Effects: | 
R3, R4, and R5 are modified by this routine. 

R6 is preserved. 


RO, R1, and R2 are used as scratch registers. RO and R1 contain the 
quadword result of EMUL that is then passed into EDIV. 


Assumptions: 
This routine assumes that all array elements Lie in the range from 0 


to 99 inclusive. (This is true if all input strings contain only legal 
decimal ergits.? The arithmetic performed by this routine will 


maintain this assumption. That is, | 
input array element LEQU 99 

times multiplier LEQU 99 | 

product LEQU 99*99 | 

plus carry LEQU 99 | 
modified product LEQU 99*100 

plus old output array element LEQU 99 

new output array element LEQU 99*101 = 9999 | 

| 


A number LEQU 9999, when divided by 100, is guaranteed to produce both 
a quotient and a remainder LEQU 99. 


a et a at a a a a a 2 ds ss 2 a 2 2 2 2 zo 
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EXTEND _STRING_MULTIPLY: 
CLRL R2 


p RO 
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; Initial carry is zero 


; Form modified product (RO LEQU 9900) 
; Add old output array elemerit 

; Remainder to output array 

; Quotient becomes carry 

; Keep going? 


ooks more complicated than it actually is. In the 
ne exits immediately. In the evens that a carr 
occurs, one additional entry in the output array w 

rare case of an output array consisting of a string of 99s will any 
ificant looping occur. 


Lt be modified, Only in 


Add final carry 

Do we overflow into next digit pair? 
Branch if carry 

Otherwise, all done 


Readjust entry and advance pointer 
Propogate or 
- and test this entry for overflow 


| 
i 
[a 
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oa » SUBTITLE VAXSDIVP = Divide Packed 


Functional Description: 


The dividend string specified by the dividend length and dividend 
address operands is divided by the divisor string specified by the 
divisor Length and divisor address operands. The quotient string 
specified y the quotient length and quotient address operands is 
replaced by the result. 


; Input Parameters: 


RO = divrlen.rw Number of digits in divisor string 
R1 = divraddr.cbd Address of divisor stcing 
R2 - divdlen.rw Number of digits in dividend string 
R5 = divdaddr.ab Address of dividend strin 
R4 = quolen.rw Number of digits in quotient string 
R5 = quoaddr.ab Address of quotient string 
Output Parameters: 
RO = 0 
R1 = Address of the byte containing the most significant digit of 
a ioe divisor string 
R35 = Address of the byte containing the most significant digit of 
R4 = 
R5 = 


the string containing the quotient 
Condition Codes: 


N <= quotient string LSS 0 
Z <= quotient string EQL 0 
: _ ge tnen overflo 


Address of the byte rag tpn Fe: most significant digit of 
Ww 


; Register Usage: 
This routine uses all of the general registers. The condition codes 
are computed at the end of the instruction as the final result is 
stares in the quotient string. R11 is used to record the condition 
codes. 

; Algorithm: 
This algorithm is the straightforward approach described in 


The Art_of Computer Programming 
Second Edition 


Volume 2 / Seminumerical Algorithms 
Donald E. Knuth 


QS awa eee 


; the dividend string 
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Notes: 


The choice of a longword array to store the auotient deserves a 
comment. In VAXSMULP, a longword array was used because its elements 
were used directly by MULP and DIVP instructions. The use of longwords 
eliminated the need to convert back and forth between longwords and 
bytes. In this roytine. the QUOTIENT_DIGIT routine returns its result 
in a register, which result can easily be stored in whatever way is 
convenient. By using longwords instead of bytes, this routine can use 
the same end processing code as MULP, a sizeable savings in code. 


- ENABLE LOCAL BLOCK 


3¢ 

3 This code path is entered if the divisor is zero. 

: Input Parameter: 

: (SP) = Return PC 

: Output Parameters: 

: O(SP) = SRMSK_FLT_DIV_T (Arithmetic trap code) 
; 4(SP) = Final state PSL 

3 8(SP) = Return PC 

: Implicit Output: 


Control passes through this code to VAXSREFLECT_TRAP. 


PD OD OF oe oe oe et ot ot ot et at et et eh ee ek ee 

PQA AAAAAAAAAA AA AAAAAAAAAAAAA AA AA AA AAA AAAAAASAAASOOSOOOOOOOOA —- x 
5O©.09 09090969 G9 0909 69 09 SII SIN NSN NAA AAA AAAAAMUMIMIUIMIUIUIUIUIN ES BEEPS EW CO 
OO NIA NEW © OO NAME WIN 9 OO NAUE WN OC OONAUE WN OWOONOUS W000 


o woo g 
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DIVIDE_BY_ZERO: 
OFFF 8F POPR #*°M<RO,R1,R2,R3,R4,R5,R6,R7,RB,RI,R10,R11> 
; Restore registers and reset SP 
7E MOVPSL <-(SP) 3; Save final PSL on stack 
04 PUSHL #SRMSK_FLT_DIV_T ; Store arithmetic trap code 
FB7D" BRW VAXSREFLECT_TRAP ; Report exception 
3; If the divisor contains more nonzero digits than the dividend, then the 
: quotient will be identically zero. Set up the stack and the registers (R4, 
3; R5, and R9) so that the exit code will be entered to produce this result. 
7E 1$: CLRL -(SP) 3; Fake a quotient digit 
59 ~=01 0 MOVL #1,R9 ; Count that digit 
FECC 1 BRW MULTIPLY _DIVIDE_EXIT ; Store the zero in the output string 
VAXSDIVP:: 
OFFF 8F 6B : 3 PUSHR #*M<RO,R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,R11> ; Save the lot 
F 91 ESTABLISH HANDLER ° ; Store address of access 
4 3 ARITH_ACCVIO ; violation handler 
é 32 ROPRAND_CHECK R4 : Insure that R4 is LEQU 31 
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49F 1 36 ROPRAND CHECK R2 3; Insure that R2 is LEQU 31 
4A7 169 MARK_ DIVP_BSBW_0 
FB56" 30 rr! ! a8 BSBW DECIMALSSTRIP_ZEROS_R2_R3 ; Strip high order zeros from R2/R3 
4AA 19 0 ROPRAND CHECK RO ; Insure that RO is LEQuU 31 
482 1701 MARK _POIN DIVP_BSBW_0 
FB4B' 30 3 \f § BSBW DECIMALSSTRIP_ZEROS_RO_R1 ; Strip high order zeros from RO/R1 
485 1704 ; Insure fhet the divisor is not zero. Because leading zeros have already 
485 1705 ; been eliminated, the divisor can gnty be zero if RO is 0 (zero Length 
485 17 $ 3; strings are ident catty zero) or 1 (R1 contains a sign digit in the low 
485 1707 ; order nibble and zero in the high order nibble). Note the an exception 
485 17 : ; will not be generated if an even length string has an oe nonzero digit 
485 1709 ; stored in its most significant nibble (including an illegal form of a zero 
re Vi ; length string. 
50 50 06 O11 EF OcBs ar; EXTZV #1,44,R0,R0 ; Convert divisor digit count to bytes 
06 12 O4BA 171 BNEQ 1 3; Skip zero divisor check unless zero 
48C 1714 MARK _POINT DIVP_0 
61 =FO 8F 3 4B8C 1715 BITB #°B11110000, TR1) ; Check for zero in ones digit 
B61 4c 1716 BEQL  DIVIDE_BY_ZERO : Generate exception if zero 
0402 1718 ; This routine chooses to do its work with a fair amount of internal storage, 
bee 1719 ; all of it allocated on the stack. The quotient is stored as it is computed, 
4C2 1720 ; in a 16-longword array. The dividend and divisor are stored as Lenguerd arrays, 
04C2 1721 ; with each array element eterna a digit pair from the orig packed 
Bee 17 ¢ ; decimal string. The numerator digits are shifted by one digit (multiplied 
4C2 1725 ; by ten) so that the quotient has its digits corregsty placed, leaving room 
04C2 1724 ; for a sign in the low order nibble of the least significant byte. A scratch 
04C2 1725 ; orrey is also allocated on the stack to accommodate intermediate results 
ber 4 $ 3 of the QUOTIENT_DIGIT routine. 
50 D6 O04C2 1728 108: INCL RO ; Include least significant digit 
58 50 7D Bees Vy 3 mova RO,R8 3; Let R8& and RO describe the divisor 
52 52 04 91 EF rte, 1731 EXTZV #1,#4,R2,R2 ; Convert dividend digit count to bytes 
§ b6 04CC 756 INCL : Include least significant digit 
7E «65 D gece Vy z MOVa R2,-(SP) ; Save dividend descriptor on stack 
56 52 50 C3 0401 1735 SUBL3 RO,R2,R6 3; Calculate main loop count 
AC IF 4) 1736 BLSSU. ‘18 ; Quotient will be zero 
56 66 e444 iy 7 INCL R6 ; One extra digit is always there 
029 1738 ; Allocate R6 longwords of zero on the stack 
50 56 DO 0409 1741 MOVL R6,RO ; Let RO be the loop counter _ 
7E 04 04DC 1248 15$: CLRL =-(§P) : Set aside another quotient digit 
FB 50 «OFS rt 4] Wt? SOBGTR RO,15$ 3; Keep going 
57 SE 00 tes 10t? MOVL SP,R7 3; Remember where this array starts 
4E4 129 ; The divisor will be stored on the stack as an array of 
4E4 1748 ; longwords. Each array element contains a number between 9 and 99 
4E4 1749 ; representing a pair of digits in the or ginal poqnes decimal string. 
4E4 17 9 3; Because the units digit is stored with the sign in packed decimal format, 
4E4 1751 ; it is necessary to shift the number as we store it. This is accomplished by 
4E4 1752 ; multiplying the number by ten. 


nm 
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; The ¢iyisor string is described by R8 (byte count) and R9 (address of most | 
: significant digit pair). 
55 8 9 C1 ADDL3 = R9,RB,RS : Point RS beyond sign digit 

2f 3 dO MOVL RB,RS 3 RS contains. the lose count 
; Put in an extra digit place for the divisor. This allows several common 

; subroutines to be used when operating on the divisor string. 

7E 04 CLRL -(SP) ; Set aside a place holder 


An orrey of longwords is allocated on the stack. R3 starts out pointing 

at the longword beyond the top of the stack. The first remainder, guaranteed 
to be zero, is ‘’stored’’ here. The rest of the digit pairs are stored safely 
below the top of the stack. 


MNEGL R8 3 ; Stack grows toward lower addresses 
MOVA (SP) n3}.SP 3; Allocate the space 
SUBL #4,SP,R ; Point R35 at next lower longword 


53 28 CE 
SE. 6E4 Dd 
53 SE 260046 COC 


MARK_POINT DIVP_R6_R7 
MOVZBL <-(R ; Get next digit pair 


GA V sd o e | 
MOVZBL DECIMAL SPACKED_T0_BINARY_TABLELRIJ.- | 
| 
| 


51 75 
51 QO00'CF41 9A 
50. 52. 251 =OA 7A 
83 52 50 00000064 8F 7 
£6 54 CF 


20S: 


7 Convert digits to binary 
EMUL #10,R1 wt ; 


: i 
br Dy ide by too" 


(R3)+ ; Divide 


EDIV R 
SOBGTR R4,20$ 


There are two cases where the final quotient (contents of R2) is zero. 

In these cases, the number of nonzero digit pairs in the divisor array is | 
smaller by one than the number of bytes containing the ortg'ags packed decimal 
oer eng- One case is a divisor string with an even number of digits. The 

second case is a divisor string with an odd number of digits but the most 
significant digit is zero (essentially a variant of the first case). The 
simplest way to handle all of these cases is to decrement R8, the divisor 
counter, if R2 is zero. Note that previous checks for a zero divisor 

prevent R8& from going to zero. 


#106 
4,2 


63 52 00 MOVL Re (R3) ; Store final quotient 
OA 1 BNEQ 25$ ; Leave well enough alone if nonzero 
56 OD INCL R6 ; One more quotient digit 
7.06 C SUBL = #4,,R7 ; Make room for it. | re - 
58 OD DECL R ; Court one less divisor ‘‘digit' 
01 = =12 BNEQ 258 


eeeeeerere BEGIN TEMP teeanneeee 


333 THE FOLLOWING HALT INSTRUCTION SHOULD BE REPLACED WITH THE CORRECT 
33 ABORT CODE. 


ze THE HALT IS SIMILAR TO THE 
fee MICROCODE CANNOT GET HERE 
323 ERRORS THAT OTHER IMPLEMENTATIONS USE. 


a a a ee ee ee a a ed a OO 
be te ee ee ee) ee > > >> >> >> > > PPP >_> >> >_>» >_>» >» ae 
BDOOOGOGDOOOS OOOO O00 OOO O00 09 09 09 09 09 09 09 09 SIN NI NINES NN NN IP AAA AAAAOAOUUIMIUIVIVVA << Oe 
DOOD NIA MEW SO OD NA UE WIN @§ 9 OD NAME WN 9 OO NOA UNE WIN = OU NOUS WN OOOO Ww 


GOOOOOOOSOOCOOSOCOCOOCOOCOOOCOOOCOOCSOSCOOOCOCOCOOOSOOOOOOOOOOOOOOCOOCOOOoOO 


SOOO OCOOOCOCOOCOM | ONUUINININNINNNININININI TO —$ DWooelovvlT TCC WOWWOMOLLrre 


a a 


B12 | 
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VAXSDECIMAL_ARITHMETIC «VAX 
v04=000 Packed =SEP=1984 00:44:34 CEMULAT.SRCJVAXARITH.MAR; 1 (7) 


halt ; This will cause an OPCDEC exception 


Ae eeeeerereee END TEMP seeeeennnee 


59 SE 00 258: MOVL SP,R9 ; RY locates low order divisor digit 
; The dividend is stored on the stack as an array of longwords. It does not 

; have its overs 08 earry snr ttee so that this storage loop is simpler. An extra 
; place is set aside in the event that it is necessary to normalize the 

; dividend and divisor before division is attempted. 


; Set aside space for U0] | 
| 


7E «4 CLRL CSF) 
52. 6746 =O! MOVAL aR YCR61,R2 ; Retrieve dividend aenerioeer 
52 62 7D MOVQ (R2) ,R2 it eee In two steps 
MARK_POINT DIVP_R6_R7 
51 83 9A 30$: MOVZBL (R3)+,R1 next decimal digit pair 
7E O000°CF41 9A MOVZBL DECIMALSPACKED_ TO BINARY. * TASLERI j,- 
492-48 P H Convert digits to binar 


From this point until the common exit path for MULP and DIVP is entered, 
no access violations that need to be backed out can occur. We do not need 


to keep the address of ARITH_ACCVIO in R10 for this stretch of code. Note 
that R10 must be reloaded before the exit code executes because the 
destination string is written and may cause access violations. 


MOVL fat ee R10 ; Retrieve size of dividend array 
MOVL SP,R11 ; R11 locates low order dividend “digit 


Allocate a scratch array on the stack the same size as the divisor array 
(which is one larger than the number of digit pairs) 


MNEGL R8,R2 ; Need a negative index 
MOVAL <-4(SP)CR2],SP ; Adjust stack pointer 


At this point, the stack and relevant pers registers contain the 


ion, N represents the number 


ed 

vi 

10 

1} 

18 

14 

iF 

f 

1§ 

5 

$5 

59 

28 y 
$3 SOBGTR R2,30$ 3; Loop through entire input string 
$) 

33 

34 

35 

36 

37 

38 

39 

40 

] 

a3 

44 

45 

rs Bey 

rk} following information. In this Seeerse 


a - Number of longwords in quotient array (M+1-N) 
- Address of pegianing of quctient arra 
- Number of dig 


LJ 
oO 
SOOOOCOCOSOOSOOOOCOCOCOCOCOCOCOOOOOSOOOSOOOSCOCSOOOOCOCOCOOCOOCOOOOOOOOOCOOCOCOCOOO =! 


00090909 09 09 09 G0 69 09 09 CD CD CD CO COCO COCO EO CO COD\WIODOOOOD OCOWOOODOODOORAOVTITO POLL LLL SSO HO 


49 ; of digit pairs in the divisor and M represents the number of digit 

29 3; pairs in the dividend. 

26 Fa a a a a a a a a a + <== SP 

os § scratch H N+1 Longwords 

54 3 Tenant anne ta pam cme mnamiet eaten + (-< R11 

22 8 dividend ' M+1 Longwords ' | 
2$ 2 terre ewe o see eco wom mene m ceca ewes + (e2< RI 
57 ; divisor H N+1 Longwords 

58 5 s+ eeeeeeeeeeeeeeeeeeeeeeeaeeeeeoe + <7 R7 

23 : quotient : M+1-N Longwords H 

61 ; H RO..R11 : 

¢ : $ ewww we ww www ween wwe ewww eo nnonn o 

Os 3 

64; 

65 ; 

66 ; 


t pairs in divisor (called N) 


| 
| 
| 
| 
| 
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VAXSDECIMAL_ARITHMETIC = VAX=11 Packed Decimal Arithmetic Instr yo i eb 9 0132348 AX/VMS Macro V04-00 Page 
v04-000 VAXSDIVP = Divide Packed -SEP=-1984 00:44:34 CEMULAT.SRCIJVAXARITH.MAR; 1 
548 1867 ; R9_ = Address of low order digits in divisor 
548 1 oe 3 R10 - Number of digit pairs in dividend (called M) 
2 : $ $ R11 = Address of low order digits in dividend 
48 1871 ° 
of F 24 1 te PUSHAL (SP) ; Store address of scratch array 
7E D Ths 187 MOVa RB,-(SP) ; Remember divisor descriptor 
7E SA 7D bee) : ie MOVQ R10,-(SP) 3; Remember dividend descriptor 
0550 1 zg ; The algorithm that guesses the quotient digit can be guaranteed to be oft. 
550 1877 ; by no more than two if the high order digit of the divisor (called V(1J) is 
550 1 of 3; at least as large as 50 (our radix divided by 2). If the high order digit 
550 1879 ; is too small, we ‘'normalize’’ the numerator and denomonator by multiplying 
beep 1380 3; them by the same number, namely 100/(vC1J+1). 
50 FC A948 01 C1 0550 13e6 ADDL3 = #1,-4(R9)CR8I,RO ; Compute VC1] + 1 
ee Boe b226 188 CMPL RO,#51 3; Compare to 50 + 1 
14 18 0559 1884 BGEQ 40 ; Skio normalization if vC£1] big enough 
53 00000064 8F 50 (C7 $228 1885 DIVL3 RO,#100,R3 ; Compute normalization factor 
54 58 7D 0563 1886 MOVQ R8,R ; Get descriptor of divisor 
QOEO 30 0566 1887 BSBW MULTIPLY_STRING ; Normalize divisor _ 
54 SA 7D 0569 1888 MOVQ R10,R ; Get descriptor of dividend 
00DA 30 056¢ 1889 BSBW  MULTIPLY_STRING : Normalize dividend 
OS6F 1891 ; We have now reached the point where we can start calculating quotient digits. 
O56F 1398 : In the following loop, R5 and R6 are loop invariants. R5 contains the number 
OS6F 1893 ; of digit pairs in the divisor. R6 always points to the lLongword beyond the 
OS6F 1894 ; most significant digit in the dividend string. R7 and R8 must be loaded on 
O56F 1895 ; each pass through because these two pointers are modified. Notice that the 
beer 1398 ; address of the divisor array is exactly what we want to store in R6. 
5A 56 7D OQS6F 1898 40S: MOVQ R6,R10 ; Let R10/R11 describe quotient and loop 
5B dD 0572 1899 PUSHL R1 ; Save quotient address for exit code 
5B 6B4A_ ODE 8228 B44 MOVAL (R11)CR103,R11 ; Store quotient digits from high end 
b278 1308 ; This rather harmless looking loop is where the work is done 
55 58 7D 0578 1904 MOva R8,R5 : Initialize count and dividend address 
59 SA 00 8278 1302 MOVL R10,R9 ; Remember the loop count in R9 
57. 10 AE =67D =O57E 1907 50S: MOVQ 16(SP) ,R7 ; Load divisor and scratch addresses 
OO1F 30 beee 1908 BSBW wad DIGIT ; Get the next quotient digit 
7B 653) «€©0 60585 «61909 MOVL R3,-(R11) ; Store it 
56 04 C2 0588 1910 SUBL : ‘‘Advance’’ dividend pointer 
FO 5A SOF bege 48 SOBGTR R10,50$ 3 «es. and go gack for more 
ate 1918 ; The quotierit digits have been stored on the stack. Eliminate the rest of the 
OS8E 1914 ; stack storage and enter the completion code that this routine shares with 
Q58E 1915 ; VAXSMULP. Note that R9 is plroady set up with the longword count used by | 
58E 1316 3; the exit code. Note also that R11 is pointing to the saved dividend descriptor 
: ; ae ; that sits on top of the saved register array. 
~~ CC ae 1919 MOVL (SP) 23° ; Reset stack pointer. . 
54 18 AB49 DE 0591 1920 MOVAL < <4#2> + - ; Skip over saved dividend descriptor 
596 1921 <4%4> >(R11)CR9],R4 |; and retrieve original R4 and R 
54 64 7D 238 13 ¢ MOVQ (R4),R4 3 eee in two steps 


41 
(17) | 
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VAX IMAL_ARITHMETIC = VAX-11 Packe —— fram Instr 16-SEP-1984 01:33:44 YAX/VMS Macro Vv04-00 Pa 4 v 
vive dts VAXSDIVP - MET met 8b i2a: 4 LCEMULAT.SR RCIVAXARITH.MAR: 1 - 9), s 
99 19 4 ; The following is @ HACK. | 
99 19 § : The sevhog used to obtain quotient digits generally leaves gerbege ee oy 
99 1927 ; in what will become the sign digit. (In fact, this is the tenths d fois of a 
599 19 8 ; decimal expansion of the remainder.) We need to make the least significant 
4 1 4 ; digit a multiple of ten. 
50 6E OA C7 0599 1931 DIVL3 #10, (SP),RO i Divide by ten. losing remainder | 
Be OA C5 HY 19 § MULL #10,0, (SP) ; Store only tens d 
FOB3=ss«331 Beh} 4 ¢ BRW MULTIPLY_DIVIDE_EXIT ; Join common exit code 
ae 1936 - DISABLE LOCAL_BLOCK 


A 
a 
A 
A 
A 
A 
A 
A 
A 
A 
A 
D 
D 
D 
D 
D 
Dd 
D 
D 
D 
D 
) 
E 
H 
s 
M 
M 
4 
M 
a 
M 
M 
M 
M 
M 
P 
Pp 
P 
P 
Pp 
Pp 
P 
Q 
S 
S 
S 
S 
S 
$ 
N 
1 
\ 
\ 
\ 
\ 
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VA TERES IAAL AAI TIME TEC = VAX-11 Packed Decimal Arithaetic Instr 16-SEP-1984 0132398 AX/VMS Macro Vv04-00 Page 43 | 
v04- QUOTIENT_DI = Get Next Digit in Quoti 5-SEP-1984 00:44:34 [CEMULAT.SRCJVAXARITH.MAR;1 (18) | 
st i . SUBTITLE QUOTIENT_DIGIT = Get Next Digit in Quotient | 
as : Functional Description: | 
a This routine divides an (N+1)-element array of longwords by an N-element | 
SA4 array, prodys ing a single quotient digit in the range of 0 to 99 
SAS inclusive. The dividend array is modified by subtracting the product 
ee of the divisor array and the quotient digit. 
5A4 The ‘‘numbers'’ that this array eperazes on multiple prqetston numbers 
SA4 in radix 100. Each digit (a number between 0 and 99) is stored in a 
5A4 longword array element with more significant digits stored at higher 
SAG addresses. The dividend string and the scratch string (also called the 
Beg product string) contain one more element than the divisor string. 
05A4 Input Parameters: 
OSA4 


RS - Number of ‘‘digits’’ (array elements) in divisor acres (preserved) 
R6 - Address of lLongword immediately following most significant 

digit of dividend string (preserved) 
R7 = Address of least significant digit in divisor string (modified) 
R8 - Address of least significant digit in product string (modified) 


ww 
> Pd 
nee 


> 
S 


Output Parameters: 


P 
1 
C 
P 
S$ 
P 
S$: 
Pp 
C 
A 


R3 = The quotient that results from dividing the dividend string 
by the divisor string. 


The final states of the three pointer registers are Listed here 
for completeness. 


mvaM— 


R6 - Address of longword immediately following most significant 
digit of dividend string 


R7 - Address of longword immediately following most significant digit 
of divisor string. This longword must always contain zero. 


R8 - Address of longword immediately following most significant 
igit of product string 


Implicit Output: 


The contents of the dividend array are modified to reflect the 
subtraction of the product string. The result of this subtraction 
could be stored elsewhere. It is a convenience to store it in the 
Sty ieene array on top of those array elements that are no longer 
needed. 


NINN NNN NA AAAAAAAAOSMIMIMIMIUIMIUIUIUIUIS SESE ELLE EWA HO 
SONA E WIN 9 OD NAU EWN CO OD NAU EW OOONOAU ES Ws" O00 


a a a a a aa a aaa a a a “ay 
S 


Side Effects: 
R7 and R8& are modified by this routine. (See implicit output List.) 
R5 and R6 are preserved. 


> 

rs 
me ek ee a ed ed = = = 3 — 0 — So = —) 2s 2s = = ss — 2 a ts ss 1 1 tt qoo 
ee PPPS EDD DE RAAAA AAR RAD $ OP PP >P>EANNNNANAANAN GT HHP P PIII HV —= 


he ty ty ty te te oy te te te a te ny ey Oe Ny ty ty ty ty ty ty ty ty ty ty ty ty ty ty ty ty ty Oy4 


000000000000 000000: 


FW CO OONOUES WN 


The contents of the divisor array are preserved. 
RO, R1, R2, and R4 are used as scratch registers. RO and R1 contain the 
| 
| 
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VAX IMAL ARITHMETIC = VAX-11 Packed Decimal Arithmetic Instr 16-SEP-1984 01:33:44 YAX/VMS Macro V04-00 P 44 | 
vive tts QUOTIENT_DIGIT = Get Next Digit in Quoti 5-SEP-1984 tree yi LEMULAT. SREIVAXARITH.MAR; 1 _ (18) 
AG 1995 ; quadword result of EMUL that is then passed into EDIV. R2 is the 
ne 1338 ; carry from one step to the next. R4 is the loop counter. | 
5a4 1998 * 
AG 1999 QUOTIENT DIGIT: 
FB AG «FC AB 00000064 BF 7A ee 2000 EMUL #100,-4(R6),-8(R6),RO ; RO <= 100 * ULj] + ULj+1) | 
50 fs A74 C SAF 01 DIVL2 -4(R7)CR5],R0 3; RO <- RO / V1) 
5 5 D 5B4 8 MOVL rh R3 3 Store quotient ‘‘digit’’ in R3 
g 1 587 0 BEQL 45§ 3 peching to do if gyet tent is zero 
00000064 8F D1 0589 4 CMPL R3,#100 : Is quotient LEQU 997 
2 1F 2 8 5 BLSSU5$ ; Branch if quotient OK 
53 00000063 8F »D0 03¢ 4 OVL #99,R3 ; Otherwise start with 99 
5C9 2008 ; We will now multiply the divisor array by the quotient digit, storing the 
B208 494 3 product in the scratch array. he n . 
2 D4 08¢9 011 5$: CLRL Re 3; Start out with a carry of zero 
ms 3 dO bece og MOVL R5,R4 3 R4 will be the loop counter 
> mm ee te back 014 10$: EMUL R3,(R7)+,R2,R0 ; Multiply next divisor digit 
88 52 50 00000064 8F 78 0503 O15 EDIV #100,R0,R2, (R8)+ : Remainder, to input array 
; Quotient becomes carr 
EF 54 FS 050¢ O17 SOBGTR R4,108 : More divisor digits? 
88 52 00 B20E O18 MOVL R2,(R8B)+ ; Store final carry 
OSE 021 ; If the product array is larger than the dividend array, then the quotient is 
OSE 0 § 3 too large. To avoid a second trip through the rather costly EMUL/EDIV loop 
OSE 023 ; and also to avoid array subtraction that produces a negative result, we will 
OSE 024 ; first compare the product and dividend arrays. If the product is smaller, we 
OSE 025 ; can safely subtract. If the product is larger, we decrease the quotient by 
B2e ° § 3 one and subtract the divisor array from the product array. | 
50 56 00 OSE 0 } 15$: MOVL R6,RO ; Point RO and R1 to high address ends 
51 58 00 OSE 0 MOVL R8,R1 3 .e. Of dividend and scratch strings | 
54 55 00 beep ? MOVL R5,R4 : Initialize the loop counter 
B268 03 ; The comparison is done from most to least significant digits 
70 71 4D OSEB 2034 208: CMPL -(R1),-(RO) ; Compare next pair of digits 
OE iF $33 035 BLSSU 30% 3; Leave loop if product is smaller 
2D 1A OSF 0 § BGTRU 50$ ; Also leave if product is larger 
F6 54 FG 2 8 i SOBGEQ R4,20$ :; More to test? 
SF 039 ; If we drop through the loop, then the dividend and product are equal. We 
SF 040 ; simply store zeros in the dividend array (the equivalent of subtraction 
SF 041 ; of equal arrays) and return. Note that RO is already pointing to the 
33 rt 3; least significant dividend array element. 
54 55 00 3 rf: MOVL R5S,R4 ; Initialize still another loop counter 
90 D4 fe £6 25$: CLRL (RO)+ ; Store another zero 
FB 54 «FG he ? SOBGEQ R4,25$ ; Keep going? 
05 O5FD rk RSB ; Return to caller 
SFE 50 
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VAX IMAL ARITHMETIC - VAX-11 Packed Decimal Arithmetic Instr 16-SEP-1984 01:33:44 YAX/VMS Macro V04-00 Page 45 | 
vive tts QUOTIENT_DIGIT = Get Next Digit in Quoti ~3Fb 1 382 0:44:34 LEMULAT. SREIVAXARITH.MAR: 1 . (18) | 
43 ; If we drop through the Loop, then the geet tone shes is stored in R3 is good. 
F ; We need to subtract the product array from the dividend array. Note that R 
FE ; and R1 need to be adjusted to point to the least significant array elements 
; : ; before the subtraction can begin. 
54 54 «OC F § 30$ MNEGL R4,R4 3; We need a negative index | 
50 6044 DE 06 MOVAL (RO)CR4].RO : Adjust dividend pointer | 
1. 61446 43 6 8 MOVAL (R1)CR4),R1 3 «ee and product pointer | 
54 55 «OD ; 25 MOVL R5.R ; R4 will count still another loop 
} 
oO HE BBE Bh pune gpe-caoe t HASTY rant FLEES 
; o end of loo no borro 
FC AO 00000064 BF C 11 8 ADDL2 #100,-4(RO) 5 Add’ borrow back to this digit : 
99 D 619 64 DECL (RO) 3 «+. and borrow from next highest digit 
EE 54 «FG rie 302 40$: SOBGEQ R4,35$ 3 Keep going? 
61E 03? : This is the exit path. R3 contains the quotient digit. The pointers to the 
Bete 308 ; various input and scratch arrays are in an indeterminate state. 
05 Bete 070 45$: RSB ; Return to caller 
061F 184 ; The first guess at the quotient digit is too large. The brute force 
O61F 2075 ; approach is to decrement the quotient by one and execute the EMUL/EDIV Loop 
O61F 2074 ; again. Note, however, that we can evaluate the modified product b 
61F 2075 ; subtracting the divisor from the initial product. Note also that, because 
6lF 878 3; the leading digit in the divisor is quarge enough’, we can only end up in 
61F 2077 ; this code path twice. (That is, the initial guess at the quotient wil 
Beir $43 3 never be off by more than two.) 
53 D7 O61F 2080 50S: DECL Ra s Ve quot tent smaller by one | 
FB S13 be 1 +t BEQL 45$ ; ALL done if zero 
06 O88 : Point R1 and R2 at the least significant digits of the scratch and product 
° bge 3 strings respectively. 
50 55 CE 06 ORG MNEGL R5,RO ; Need a negative index 
51 FC A840 DE 0626 2087 MOVAL rrtiy fafa R1 3; Scratch array contains N+1 elements 
52. 6740 43 628 +43 MOVAL (R7)CROJ,R2 ; Product array contains N elements 
54 55 OD > F 434 MOVL R5, ; R4 will count still another loop 
81 Se 6 44 60$ SUBL2 (R2)+,(R1)+ ; Subtract next digits 
A 1 6 35 BGEQ ; Skip to end of loop if no borrow 
FC Al 00000064 8F CO 06 09 ADDL2 #100,-4(R1) ; Add borrow back to this digit ‘ 
1 D7? 063F Boe DECL (R1) 3 «+. and borrow from next highest digit 
EE 54 FG 2) 95 70$ SOBGEQ R4,60$ 3; Keep going? 
51 04 CO 0644 099 ADDL2 #4,R1 ; Point R1 at most significant digit 
99 11 0647 98 BRB 15$ ; Make another comparison 
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VaERES SSEISEE Ob:dEt wo tf 


ply a String by a SEP- 0:44: EMULAT.SRCJVAXARITH.MAR; 1 
. SUBTITLE MULTIPLY_STRING = Multiply a String by a Number 
Functional Description: 


Pa 
ST 


‘io 
o 
zo 

c- 
-23 
re 


This routine mg he an array of numbers (each array element LEQU 
9) by a number (also LEQU 99). Each array element in the input array 
is replaced with the modified product, with the carry propogated to 
the next array element. 
Input Parameters: 
R3 - Multiplier 
R4 = Input array size 
RS - Input array address 
Output Parameters: 
None 
Implicit Output: 
The input array elements are altered. 
Side Effects: 
R4 and RS are modified by this routine. 


R3 is preserved. 


PAA AAAPAASASAOA SOSA OSOSASSOSS r— 
De ee ee cee 
WDOOOOOOVOOOOOOVOOOOOOOOOOOO | 


RO, R1, and R2 are used as scratch registers. RO and R1 contain the 
quadword result of EMUL that is then passed into EDIV. RZ is the 
carry from one step to the next. 


Assumptions: 


This routine assumes that all array elements Lie in the range from 0 
to 99 inclusive. (This is true if all input strings contain only legal 
decimal digits.) The arithmetic performed by this routine will 
maintain this assumption. The details of this argument can be found in 
the routine header for EXTENDED _MULTIPLY_STRING. This routine performs 
less work so that those arguments also apply here. 


MULTIPLY_STRING: 
CLRL R2 


RO ODNA NE WIN 9 OD NAME WIN 9 OO NAME WIN 0 OD NAME WN OODNAUMUESWNO Oa 


UAE REESE EES WWW NWI NPN  BOODOOOOOOOCO Ze 


a a ad et ot a 8 4 ss ss 5 ss st a tt —z 


52 D4 oH ; Initial carry is zero 
50. 652) 6650S 537A (064 108: EMUL R3,(R5),R2,RO ; Form modified product (RO LEQU 9900) 
85 52 50 00000064 8F 7B $30 EDIV #160,R0.R2.(R5)+ 3; Remainder to input array 
65 3; Quotient becomes carry 
EF 54 «FS 22? SOBGTR R4,10$ ; Keep going? 
65 52 0 65C MOVL R2,(R5) ; Store final carry 
0 65F RSB 


a= 
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VAX IMAL ARITHMETIC Decimal Arithmetic Instr 16-SEP-1984 01:33:44 YAX/VMS M 
vou~000 Bn 3Ep=198e bosee $e RCI 


ro V04-00 
SEP- EMULAT.SRCJVAXARITH.MAR; 1 


é ] : j . SUBTITLE DECIMAL _ROPRAND 
5 3 Functional Description: 

4 158 ; This routine receives control when a digit count larger than 31 
6 1 ; is detected. The exception is architecturally defined as an 
6 160 ; abort so there is no need to store intermediate state. ALl of the 
66 161 ; routines in this module save all registers RO through R11 before 
6 196 3 ot begs | the digit check. These registers must be restored 

re 197 ; efore control is passed to VAXS$ROPRAND. 

29 193 ; Input Parameters: 
6 189 : 00(SP) = Saved RO 

$60 9169 | 
6 19 : 44(SP) = Saved R11 

re 17 3 48(SP) = Return PC from VAX$xxxxxx routine 

des 178 : Output Parameters: 

0660 2175: OO(SP) - Offset in packed register array to delta PC byte 

Bee8 176 3 04(SP) = Return PC from VAX$xxxxxx routine 

dee 178 + Implicit Output: 

bee 139 ; This routine passes control to VAXS$ROPRAND where further 

ret 183 3 exception processing takes place. 

660 2183 © 

asso 184 ASSUME ADDP6_B_DELTA_PC EQ ADDP4_B_DELTA_PC 

0660 2185 ASSUME SUBP4"B-DELTA_PC EQ ADDP4~B-DELTA_PC 

660 21 6 ASSUME SUBP6"B-DELTA-PC EQ ADDP4~B-DELTA_PC 

66 1 ASSUME MULP_B_BELTA_PC EQ ADDP4"~B-DELTA_PC 

se) : 3 ASSUME DIVP_B_DELTA_PC EQ ADDP4_B_DELTA_PC 

3088 190 DECIMAL_ROPRAND: 

OFFF 8F BA 660 191 POPR #*A<RO,R1,R2,R3,R4,R5,R6,R7,RB,RI,R1I0,R11> 
03 bp 664 198 PUSHL #ADDP4’B_DELTA_PC : Store offset to delta PC byte 
F997" 31 0666 219 BRWw VAX$SROPRAND ; Pass control along 
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VAX IMAL ARITHMETIC = VAX-11 Packed Decimal Arithmetic Instr 16-SEP-1984 01:33:44 VAX/VMS Macro v04-00 Page 48 Vv 
vOee000 ARITH_ACCVIO = Reflect an Access Violati 5-SEP-1984 $b iaz:4e YEMULAT. SREIVAXARITH.MAR: 1 s (21) vl 
é 193 na - SUBTITLE ARITH_ACCVIO = Reflect an Access Violation 
199 : Functional Description: 
66 35 : This routine receives control when an access violation occurs while 
669 ; executing within the emulator routines for ADDP4, ADDP6, SUBP4, SUBP6, 
424 ; MULP, or DIVP. 
$3 : : The routine header for ASHP_ACCVIO in module VAXSASHP contains a 
69 4; detailed description of access violation handling for the decimal 
9 5 ; string instructions. 
$63 : : Input Parameters: 
94 09 ; See routine ASHP_ACCVIO in module VAXSASHP 
re 1} : Output Parameters: 
289 18 ; See routine ASHP_ACCVIO in module VAXSASHP 
669 2215 ° 
669 1 ARITH_ACCVIO: 
52 D4 0669 1 CLRL oa : Initialize the counter 
F991 CF OF 0668 18 PUSHAB MODULE_BASE ; Store base address of this module 
O6D3'CF 9F O066F 1 PUSHAB MODULE END 3; Store module end address 
F9BA’ 30 0673 0 BSBW DEC IMACSBOUNDS_ CHECK ; Check if PC is inside the module 
5E 04 3} 978 1 ADDL #4,SP 3; Discard end address 
51 BE OC 679 : SUBL2 (SP)+,R1 : Get PC relative to this base 
O0OO'CF42 51 B1 O067C 4 108: CMPW R1,PC_TABLE_BASECR2) ; Is this the right PC? 
07 13 0682 5 BEQL ; Exit loop if true 
F452 28 Fe tt § AOBLSS #TABLE_SIZE,R2,10$ ; Do the entire table 
ons 8 ; If we drop through the dispatching based on PC, then the exception is not 
3 : 3 3 one that we want to back up. We simply reflect the exception to the user. 
OF BA és 1 208: POPR #*M<RO,R1,R2,R3> 3; Restore saved registers 
05 ° A ¢ RSB 3; Return to exception dispatcher 
re 4 ; The exception PC matched one of the entries in our PC table. R2 contains 
oe 5 ; the index into both the PC table and the handler table. R1 has served 
§ ; its purpose and can be used as a scratch register. 
51 Q000°CF42 75 68 8 30$: MOVZWL HANDLER_TABLE BASECRZI A ; Get the offset to the handler 
F96A CF41— 1 at i, JMP MODULE _BASECRT 3; Pass control to the handler 
D696 41 ; In all of the instruction-specific routines, the state of the stack 
696 g 3; will be shown as it was when the exception occurred. ALl offsets will 
0696 43 ; be pictured relative to RO. 
| 
} 
| 


——.. ee Sa 1 


K 12 
VAX IMAL ARITHMETIC = VAX-11 Packed Decimal Arithmetic Instr 16-SEP-1984 01:33:44 YAX/VMS Macro V04-00 P 49 
vive tt Access Viotstion Mandl ine for AGOPn nnd” SSep-ioRe Ooseccde PERU AE. caeScavaRisneman:1 P29 (33) 


45 cs . SUBTITLE Access Violation Handling for ADDPx and SUBPx 
; Functional Description: 


| 
The only difference anene the various entry points is the number of | 
Lon words y he stack. RO is advanced beyene these longwords to point | 

sto | 


696 
69 6 
$96 $52 
$9 : 
69 : he | saved registers. These registers are then restored, 
24 : effectively backing the routine up to its initial state. 
2 2¢ ; Input Parameters: 
$38 2$ ; RO - Address of top of stack when access violation occurred 
44 28 ; See specific entry points for details 
Boge ¢ : Output Parameters: 
ooe 6¢ : See input parameter List for VAX$DECIMAL_ACCVIO in module VAXSASHP 
ooe $504 
0696 66 ; ADD_SUB_BSBW_24 
0696 68 : An access violation occurred in one of the subroutines ADD_PACKED_BYTE, 
0696 69 ; SUB_PACKED_BYTE, or STORE_RESULT. In addition to the six longwords of work 
0696 70 ; space, this routine has an additional longword, the return PC, on the 
beoe 4 3 stack. 
0696 73; 00(RO) = Return PC in mainline VAX$xxxxxx routine 
0696 74; 04(RO) - Address of a byte of destination string 
0696 75 3 08(RO) - First longword of scratch space 
0696 76 ; etc. 
Bb Bere 
0696 ig ADD_SUB_BSBW_24: 
50 04 CO be98 + ADDL #4,R0 ; Skip over return PC and drop into ... 
0699 3+ 
9944 5 3; ADD_SUB_24 
0699 5 ; There are five longwords of workspace and a saved string address on the stack 
+944 : for this entry point. 


) - First longword of scratch space 


- Fifth longword of scratch space 
d RO 


0(RO) 

4(SP) Save 
8(SP) = Saved R1 
etc. 


| 
‘oo i 
j 


699 3 : sett - Address of ote byte of destination string 
3 4(RO 


9 
699 38 ADD_SUB_24: , 
18 69 699 9 ADD #24,R0 ; Discard scratch space on stack — 
F961" 31 te 00 BRW VAX$DECIMAL_ACCVIO ; Join common code to restore registers 
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F 

F 
e9F 
69F 
69F 
69F 
69F 
69F 
69F 
oor 
Deer 
69F 

50 04 C0 


069F 
F95B" 31 O6A2 23 


oatty point has an additional longword, the return PC, on the stack on top 
of the saved register array. 

O(RO) = Return PC in mainline VAX$xxxxxx routine 

4(RO) = Saved RO 

oor - Saved R1 

etc. 


: ADD _5UB_BSBW_0 
; An access violation occurred in one of the subroutine STRIP_ZEROS. This 
1 
1 
1 
1 
1 
1 
1 
1 


0 
1 
§ 
4 
5 
6 
? 


ADD_SUB_BSBW_0: 
ADDL #4,R0 ; Skip over return PC and ... 
BRW VAXSDECIMAL_ACCVIO ; Join common code to restore registers 


. 4 


Lee 
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- SUBTITLE Access Violation Handling for MULP and DIVP 


19 
Q ;+ 
oA 1 : Functional Description: 
6A : : The only difference enone the various entry points is the number of 
6A 4; lLongwords on the stack. RO is advanced peyone these longwords to point 
6A > 3 to the list of saved registers. These registers are then restored, 
oA § i effectively backing the routine up to its initial state. 
* 8 : Input Parameters: 
oA 9 ; RO - Address of top of stack when access violation occurred 
Deh § : See specific entry points for details 
Bene : ; Output Parameters: 
pone § : See input parameter List for VAX$DECIMAL_ACCVIO in module VAXSASHP 
BA Ss 
O6A5 ¢340 > MULP_R8 
06A5 ri} ; An access violation occurred while MULP was accessing one of its two source 
06A5 45 ; strings. In this particular case, MULP was eter tas the longer of the two 
06A5 44 ; input strings in a longword array on the top_of the stack. There is an 
06A5 45 ; array of R8& longwords on top of an array of 32 lLongwords on top of the 
BoA? 4) 3 saved register array. 
Bons rt: : R8 - Number of longwords on top of the 32-lLongword array 
06A5 2350 * 
06A5 331 ENABLE LOCAL_BLOCK 
O6as 5385 MULP_RB 
50 6048 ODE O06A5 54 ~""MOVAL (RO) CRB, RO ; Discard input array storage 
04 11 4 56 BRB 10$ ; Might as well share a Little code 
06AB 57 ;+ 
aan 335 : MULP_AT_SP 
deAB 60 : An access violation occurred while MULP was accessing one of its two source 
6AB 61 ; strings. In this case, the access violation occurred in the middle of the 
6AB 6¢ 3 grand multiply loop as a digit pair was being retrieved from the shorter of 
6AB 65 ; the two input strings. The address of the start of the 32-longword array 
7 rt: 3 was itself stored on top of the stack for convenience. 
6AB 66 : 90¢R0) - Saved byte count of longer input string 
6AB 67 ; 4(RO) = Saved aserons of longer input string 
7 $3 3 08(RO) - Address of 32-longword array farther down the stack 
6AB 0° 
6AB 2371 MULP_AT_SP: 
50 08 ap DO 06AB % MOVL 8(RO) ,RO ; Locate start of 32-longword array 
50 0088 C 35 6AF 73 108: MOVAB <<4#30> + <4#2>>(RO),RO : Throw that away, too ; 
F949" oBs ee BRW VAXSDEC IMAL_ACCVIO ; Join common code to restore registers 


a a iia a 7 
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Saved R1 


where zz = 4 * (R9 = 1) 


Access Violation Handling for MULP and D 5-SEP-1984 EMULAT.SRCJVAXARITH.MAR; 1 
$31 6 - DISABLE LOCAL BLOCK 
78 ;+ 
rtf 8 3; MULP_DIVP_R9 
687 1 ; An access violation occurred while the final result was being stored in the 
687 ¢ ; result string. In this common exit code path, R9 counts the number of | 
687 ; longwords on the stack. In all cases where an tay violation can occur, a | 
687 4 ; longword has been removed from the stack but R9 has not yet been | 
687 5 ; decremented to reflect this. The conceptual instruction sequence that 
687 § 3 resets the stack pointer (really RO) to point to the start of the saved 
687 3; register array is 
687 8 ; 
687 9; DECL Rg 
687 0; MOVAL (RO)CR9I 
687 91; 
oe? ; ; A single instruction accomplishes this. 
687 94 ; R9 = One more than the number of longwords on the stack on top 
af 32 3 of the saved register array. 
oat 344 3 0O(RO) - First Longword of scratch storage remaining on the stack 
687 2399 ; . 
687 2400 ; zz-4(RO) = Last longword of scratch stereer : 
687 2401 ; zz+0(RO) = Saved count of dividend or multiplier etring 
687 2402 ; zz+4(RO) = Saved address of dividend or multiplier string 
687 2403 ; zz+8(RO) = Saved 
687 3 2(RO) = 
687 3 
687 3 
687 ; 
687 : 
687 

687 MULP_DIVP_R9: 
50 04 A049 DE 0687 MOVAL 4(RO)CR9],RO ; Discard scratch storage on stack 

F941" 31 ree) BRW VAXSDECIMAL_ACCVIO ; Join common code to restore registers 

68 + 

68 MULP_DIVP_8 

68 

68 


SOOOOCCOCOCOCOOCOSCOOOCOOCOOOOOOOOOOCOoO 


ee De Ds De | 


An access violation occurred in the common exit path after the scratch array 
had been removed from the stack but before the saved descriptor for the 
multiplier string was discarded. 


SSSSS 

ee es Bes Bee Be | 
OOP EOC E FEELS RAAT Tt = 
ROOD UNSW) "OVO OWNOU SW "0 0On 


GAINIWNRORoRoRoRoNORoNofons 3 a 20 


O(RO) = Saved count of dividend or multiplier string 
4(RO) = Saved address of dividend or multiplier string 
6BF 8(RO) = Saved RO 
6BF 12(RO) = Saved R1 
6BF etc. 
O68F - 
O6BF 
6BF MULP_DIVP_8 i ren ; . 
50 98 0 68F ABD #8,RO ; Discard multiplier string descriptor 
F93B" 31 ecs BRW VAXSDECIMAL _ACCVIO ; Join common code to restore registers ' 
6C i+ 
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PERLE LL SRL RRR RRP PP PPP PPP PP PPP PPP PP PPE 


009090909 00 09 09 09 09 09 SI NINN SN NN SND PA AAA AAA A MTT PPP EEE 


POPIPOPIPOPIPININININPINININININPINIPYINIPININIPONY 


0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
: 


PREPARE DPAEDPPAEDPRDPPA AAA AAA AAPA AA AAA AA AAA AA AA AAA AA AAA AAAAAAAAOA OAS 
PAAAADAAAAOVOAOAAOAAOAAAAAAAAAAAAAAAAAAABAAAAHAHAAOHAAAOOAHOAAOOMOOO 
DBDWOWWWOWVOOWOOWVOOOOOW WWW WW W 00000000 00 0000 09 09 09 09 09 09 CO COCO COMI IMI IIIT 
DOONAN EWR CO OO NAUE WIND 0 OD NAW EWI 3 O OO NA UNE WW OOO NOAUE WR) OOONOAUEwW 


PPP PS SS SS SS SS SS SS 


POPIPONININPYONPINPIPNONIPONIPONI 


An access violation occurred in one of the subroutine STRIP_ZEROS. This 
eatty point has an additional longword, the return PC, on the stack on top 
of the saved register array. 
O(RO) = Return PC in mainline VAXS$MULP or VAX$DIVP routine 
4(RO) = Saved RO 
re - Saved R1 
etc. 


ULP_BSBW_0: 
IVP~BSBW~0: 
ABDL #4,R0 ; Skip over return PC and drop into ... 


DT Pe Se Se Se Se Se Ge Gee. SH Be ee 


+ 


DIVP_O 
MULP_DIVP_O 
There was nothing allocated on the stack other than the saved register 


array when the access violation occurred. We merely pass control to common 
code to restore the registers. 


0O(RO) = Saved RO 
ie - Saved R1 
etc. 


cue 
cr-< 
‘hs hes 
T=) 


P_0: 
BRw VAX$SDECIMAL_ACCVIO ; Join common code to restore registers 


DIVP_R6_R7 


An access violation occurred while one of the two input strings was being 
converted to an array of longwords_on the stack. The state of the stack 
is rather complicated but R6 and R7 contain enough information to allow 
the rest of the stack contents to be ignored. 


R6 = Count of Longwords in quotient array on stack 
R? - Address of quotient array on stack 


00(RO) - First Longword of quotient array 


22z-4(RO) = Last longword of scratch storage 
22+0(RO) = rigs count of dividend string 
22+4(RO) - Address of dividend string 
zz+8(RO) = Saved RO 

22+12(RO) - Saved R1 

etc. 


where 22 = 4 * R6 


Be Se Oe Oe Se Se Ge Se Be Be Be Ge Se Ge Se Ge Ge Ge Ge Ge Ge Be Se 
‘ + 


44 \VAX/VMS Macro V04-00 Page 
g for MULP and D 5-SEP=-1984 00:44:34 (CEMULAT.SRCJVAXARITH.MAR; 1 
MULP_BSBW_ 
DIVP_BSRW_ 


i 
| 
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6cB 
6CB a3! DIVP_R6_R7: 
50 08 A746 BF of 13% “MOVAL 8(R7)CR61,R0 ; Discard eoreens on stack 
F92D" rt rh BRW VAXSDECIMAL_ACCVIO ; Join common code to restore registers 
6D 495 END_MARK_POINT 
06D 496 
06D 497 «END 
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sacPhecs = ° 328 VAXSDECIMAL_EXIT eeeeeeee § X 0 
.»»ROPRAND... = 499 R 02 VAXSDIVP 90000488 RG $ 
ADDP4_B_DELTA_PC = 000 H VAXSMULP_ __ 0000287 RG 
ADDP6_B DELTA_PC = 0000000 VAXSREFLECT_TRAP areneree = X 00 
ADD_PACRED DA R 0 VAXSROPRAND eeeeenee =X 00 
ADD-PACKED_BYTE_R6_R7 3 5 R 8 VAXS$SUBP4 99000082 RG 02 
ADD- PACKED BYTE STRING QO0015F R VAX$SUBP6 0000000 RG 02 
ADD" SUBTRACT _EXTT 0001 3 R 0 

“SUB_24 0699 R 0 
ADD-SUB_BSBW_ O69F R 0 
ASS aus-Y-PEROcR tt an 

= 

ARITH_ACCVIO ~ $o000869 B f3 
DECIMALSBINARY_TO_PACKED_TABLE teeeeeee x (0 
DECIMALSBOUNDS_ CHECK eeeeeeee =X 00 
DECIMALSPACKED_TO BINARY TABLE teeeeeee x (0 
DECIMALSSTRIP_ZEROS_RO_RT seeeeeee =X 00 
DECIMALSSTRIP~ZEROS_R2_R3 eeeeeere =X 00 
DECIMAL_ROPRAND 00000660 R 02 
DIVIDE_By_ZERO 00000478 R 02 
DIVP_ 444943 4 R 02 
DIVP_BSBW_O 000006C5 R 02 
DivP_B DECTA_PC = 00000003 
DIVP-RE_R 000006CB R 02 
EXTEND STRING_MULTIPLY 0000044A R 02 
HANDLER_ TABLE _ 00000000 R 04 

ULE _BASE = 00000000 R 02 
MODULE _END = 000006D3 R 02 
MULP_AT_SP 000006AB R 02 
MULP~BSBW 000006C5 R 02 
MULP~B_DECTA_PC = 00000003 
MULP~DTVP_ 000006C8 R 02 
MULP_DIVP_8 QOOOO06BF R 02 
MULP_DIVP_R9 00000687 R 02 
MULP_R seeie rtd fe 02 
MULTIPLY_DIVIDE_EXIT 00000557 R 02 

TIPLY_STRING 00000649 R O¢ 
PC_TABLE_BASE 00000000 R 0 
PSCSM_N = 00000008 
PSLSM_V = 00000002 
PSLSM_ 5 4 4 4 
PSL$V_CM = 0000001F 
ratty = $eeetee 
QUOTIENT _DIGIT BO00SAG @ 02 
SRMSK_FLT _DIV_T = at ye: 
STORE_RESOLT Boe 49 R 02 
SUBP4"B_DELTA_PC = 00000 9 
UBP6-B" DELTA _PC = 000 
SUBTRACT PACKED 00001 3 R 9 
SUB_PACKED_BYTE_R6_R7 4 R 
SUB” PACKED_BYTE_ STRING 0 1D R 0 
TABLE SIZE = ° 8 
VAXSADDP4 B RG 0 
VAXSADDP6 9 RG 8 
VAXSADD_PACKED_BYTE_R6_R7 165 RG 
VAXSDECTMAL_ACTVIO eeeerere =X 0 


VAXSDECIMAL _ARITHMETIC 
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PSECT name Allocation PSECT No. Attributes 


- ABS . 0000000 ( -) 0 ¢ 0.) NOPIC  USR ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
SABSS 000000 ( -) ) NOPIC USR CON ABS LCL NOSHR’~ EXE RD WRT NOVEC BYTE 
VAXSCODE 0 Boee ( 1747.) ( ¢-} PIC USR CON REL tL HR EXE RD NOWRT NOVEC LONG 
C_ TABLE 4 005 ( 68) 8 ( ) PIC USR CON REL LCL SHR NOEXE RD WNOWRT NOVEC BYT 
HANDLER_TABLE 0000005 ( 0.) 4¢ 4.) PIC USR REL LCL SHR NOEXE RD NOWRT NOVEC BYTE 

TT 

' ; Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 10 00:99 :00.06 00:00:00.99 
Command processing 71 00:00:00.5 00:00: 93 -¢¢ 
Pass 1 208 00:00:07.77 00:00:22.36 
Synbol. table sort 0 00:00:00.35 00:00:01.58 
Pass 2 392 00:00:04.76 00:00:13.45 
Symbol table output 0 bs Ba 00:00: +0) 68 
Psect synopsis output 0 00:00:00.0 00:00:00.0 
Cross-reference output 0 aes Se 00:00:00.00 
Assembler run totals 681 00:00:13.58 00:00:42.27 


The pereins set Limit was 1650 pages. 

50325 bytes (99 pages) of virtual memory were used to buffer the intermediate code. 

There were 20 pages of symbol table goece ai Secotes to hold 182 non-local and 113 local symbols. 
2497 source Lines were read in Pass 1 "aie he object records in Pass 2. 

23 pages of virtual memory were used to define 21 macros. 


woe et wee eom ese ese eee ee eee et 


Macro Library name Macros defined 


-$255$DUA28 : CEMULAT .OBJ J VAXMACROS.MLB; 1 12 
$255$DUA28: SYSLIBISTARLET. MLB;2 6 
TOTALS (all Libraries) 18 


318 GETS were required to define 18 macros. 
There were no errors, warnings or information messages. 
MACRO/LIS=LIS$:VAXARI TH/OBJ=0BJ$:VAXARITH MSRC$:VAXARI TH/UPDATE=(ENHS$: VAXARITH) +LIB$: VAXMACROS/LIB 
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